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PREFACE 

The simpler observations in X-ray crystallography can 
be treated on the idea that X-rays are reflected from 
layers of atoms in the crystal, and that the beams 
reflected by successive parallel layers interfere with 
one another. This was sufficient, for example, for 
the interpretation of the early ionization-spectrometer 
measurements, and of the Laue photographs of simple 
substances. When more complex crystals were 
studied some means of exhibiting neatly the informa- 
tion given by the large numbers of spots on moving- 
crystal and moving-crystal moving-film photographs 
were necessary, and the device known as the reciprocal 
lattice was found convenient for this purpose. This 
device is indeed almost a necessity for the interpreta- 
tion of such photographs. The reciprocal lattice is 
not generally familiar, and many books on crystal- 
lography stop short at the point where its use becomes 
desirable. This point is roughly the division between 
one-dimensional and three-dimensional problems in 
X-ray optics. Some of the latter problems are treated 
here : the construction and properties of the reciprocal 
lattice, its relation to the commoner types of X-ray 
diffraction photograph, the representation in reciprocal 
space of crystals that are small, distorted or otherwise 
imperfect, and the effects of the imperfections on 
-^•-ray photographs. The chapter on the effects of 
thermal agitation on X-ray diffraction is included with 
s °rne misgivings, since the writer has no first-hand 
a cquaintance with it. One subject that would natur- 
a % form part of X-ray Optics is omitted : the use of 
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Fourier series in the investigation of crystal structures. 
The omission is partly because of the considerable 
space that would be required to deal with it usefully, 
and partly because, in the no-man's-land of Crystal- 
lography, structure determination is definitely closer to 
Chemistry than to Physics. 

It is hoped that a satisfactory compromise has been 
struck between the diffuseness of a non-mathematical 
treatment and the uninviting symbolism of some of 
the original papers. The necessary mathematics^ is 
introduced gradually, the earlier chapters requiring 
little more than elementary trigonometry. The mono- 
graphs X-rays and X-ray Crystallography in this series 
have been kept in mind, and it will be found that there 
is practically no overlapping in treatment and very 

little of topic. 

A. J. C. W. 

UNIVERSITY COLLEGE 
CARDIFF 

September 1946 
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CHAPTER I 

INTRODUCTION 

i. Ideally, a crystal * is an arrangement of atoms that 
repeats itself indefinitely in three dimensions. The repeat- 
ing group of atoms occupies a small parallelepiped, called 
the unit cell. The edges of the unit cell are the crystal 
axes, and are generally represented by the letters a, b, c. 
The parallelepiped is not necessarily rectangular, and the 
angles between the pairs of axes arc denoted cr.,(i,y (figure i). 




Fig. i. — Axes and angles of a unit cell. 

The complete crystal is made by the regular packing to- 
gether of these units, and the network formed by the edges 
of the unit cells is called the crystal lattice. In most 
crystals the atoms within the unit cell are arranged accord- 
ing to certain geometrical rules of symmetry that govern 
the positions of equivalent atoms. It can be shown that 
there are exactly two hundred and thirty sets of rules 
Possible. They are known as space groups, and are of 
fundamental importance in crystal structure analysis. A 
f eally perfect crystal would be of infinite size, as the 
Periodicity necessarily stops at the crystal boundaries, but 

Fuller descriptions and more precise definitions of the 
ec hnical terms introduced in this paragraph will be found, for 
^mplc, in The Crystalline State, by W. L. Bragg (Bell, London, 
9 33), or in X-Ray Crystallography, by R. W. James. 
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for most purposes a crystal of more than a thousandth of 
a centimetre in linear dimensions can be considered perfect. 
The main interest of this book is in the diffraction effects 
produced in X-ray photographs of a crystal in which some- 
thing has gone wrong — its size too small, or the lattice- 
network curved, or the space-group rules only approximately 
obeyed. 

2. It is convenient to begin by recalling some results 
that can be obtained without the use of the reciprocal 
lattice. The first of these is Bragg's law. In deriving 
this law it is postulated that a crystal consists of parallel 
planes of atoms whose arrangement repeats at regular 




Fig. 2. — Derivation of Bragg's law and estimation of the angular 
range of reflection. 

intervals d (figure 2). Such atomic planes arc very evident 
when a model of some simple crystal is examined, and may 
be picked out in many different ways. These planes are 
each supposed to reflect a certain small fraction of the 
incident beam of X-rays, so that the angle between the 
incident rays and the planes is equal to the angle between 
the reflected rays and the planes, and so that the incident 
and reflected rays lie in the same plane. For an arbitrary 
angle the reflected intensity is small, as there is then no 
special phase relation between the reflections from succes- 
sive planes, and the reflected rays largely cancel one an- 
other. For certain angles, however, the extra distance 
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travelled by the ray A'B'C will be an integral number of 
wave-lengths greater than that travelled by the ray ABC, 
and the reflected rays reinforce one another. It is easy to 
see that, if this condition is fulfilled for the pair of rays 
ABC and A'B'C, it is fulfilled for all other rays, and a 
strong reflected beam is produced. If BD and BE are the 
perpendiculars from B on A'B' and C'B', the extra distance 
travelled is DB' + B'E, and reinforcement occurs if this 
distance is n?., where n is any integer and X is the wave- 
length of the monochromatic incident X-rays. The lengths 
of DB' and B'E are each d sin 0, so the beam is strongly 
reflected when 

2d sin = n?., . . (i) 

which is Bragg's law. The deviation, the angle between 
the directions of the incident and reflected rays, is twice 
the Bragg angle 6. The reflection with n = i is called the 
first-order reflection, that with n = 2 is the second order, 
and so on. 

3. A simple extension of this argument, suggested to the 
writer by Dr A. R. Stokes, leads to an estimate of the 
angle /3 over which the intensity of reflection is appreciable. 
At the angle given by equation (1) all the planes are reflect- 
ing so that their contributions add, but the contribution of 
the last plane is actually mn wave-lengths out of phase with 
that of the first. If it were mn + 1 or mn - 1 wave-lengths 
out of phase the intensity of reflection would be zero, since 
the contribution of the top half of the crystal would be 
cancelled by the contribution of the bottom half. As a 
rough estimate the range of over which the intensity of 
reflection is appreciable may be taken as half the difference 
between the angles giving these phase relationships. If 

2£ sin 0' = {mn + i)l, 
2£ sin 6" = (mn - i)A, 
e(sin0' -sinO") - X,. 

6' +6" . 0' -0" , 
2£ cos sin - = /., 
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or, since 0' and 0" differ very little from 0, 
0' -0" =X/ecosQ, 

where s = md is the thickness of the crystal. If e is small 
(less than about io~ 4 cm.), ($' - 0") is quite appreciable, 
and the crystal reflects over a small range on either side 
of the angle given by equation (i). The effect is analogous 
to the imperfect resolution given by diffraction gratings of 
few lines. Since the actual deviation of the rays is 20, 
the range of deviation over which the intensity is appreciable 
is (0' - 0"), or 

P~A/eco&0. . . . (2) 

This argument gives only the order of magnitude of /5, 
but the proportionality factor is not far from unity, and 
it is convenient to define the ' apparent particle size ' of 
a crystal by 

E^X/ficosO. ... (3) 

The intensity of reflection does not, of course, remain 
constant over an angular range /? and then drop suddenly 
to zero. A more precise definition of /5 in terms of the 
measured variation of intensity of reflection with angle 
is given in chapter II, and its variation with the shape 
and orientation of the crystal particle is discussed in 
chapter IV. 

4. If the crystal is larger than io -4 cm., reflection may 
still take place over a considerable range if the crystal is 
otherwise imperfect. For example, in a distorted crystal 
the interplanar spacing d will not be constant throughout, 
and different parts of the crystal will reflect at different 
angles. The broadening of the lines on powder photo- 
graphs of an aggregate of such distorted crystals may be 
estimated from equation (1). If the maximum strain in 
the crystal is e, the spacing will vary from d(i + e) to 
d{\ - e), and will vary over a range numerically equal 
to 2ed(d0/dd), or ze tan 0. Since the deviation of the rays 
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is 20, the range of angle over which the reflection is appreci- 
able will be twice this, or 

£~4«tan0. . . . (4) 

In cases of broadening from this cause it is convenient to 
define the ' apparent strain ' by 

rj = /S cot 0. . . . (5) 

Equation 4 suggests that the apparent strain is about four 
times the maximum strain. The relation between actual 
and apparent strain is treated in chapter VIII, but the 
problem is not in as satisfactory a state as that of the 
broadening due to small crystal size. 

5. The other forms of crystal imperfection do not lend 
themselves to elementary treatment, and it seems best to 
postpone discussion of them to later chapters. 






CHAPTER II 

X-RAY DIFFRACTION AND THE RECIPROCAL 
LATTICE 

i . The derivation of the Bragg law given above is simple, 
but it has certain weaknesses. Aesthetically it is displeas- 
ing, because it is a mixture of geometrical and physical 
optics. The interaction of the X-rays with a single plane 
of atoms is treated by geometrical optics, as simple reflec- 
tion, while the interaction with successive planes is treated 
by physical optics, as interference. It is by no means 
obvious why a single plane of atoms ought to reflect X-rays 
as a mirror reflects light ; we shall see shortly that in fact 
it would not. Also, in complex structures, where there 
are atoms in ' general ' positions in the unit cell (that is, 
atoms whose positions are not determined completely by 
symmetry) the successive planes of atoms postulated are no 
longer evident, and the treatment loses its plausibility. We 
must therefore make a new start on the problem of the 
interaction of X-rays and crystals, using physical optics 
throughout, and treating the whole crystal impartially, with 
no bias toward planes and interplanar spaces. This will 
lead naturally to the important concepts of reciprocal space 
and the reciprocal lattice, which form the basis of modern 
X-ray crystallography. 

2. Suppose first that each unit cell contains a single 
scattering centre. This may be imagined as a single atom, 
or as the centre of gravity of the actual distribution of 
electrons in the unit cell. The fact that all the electrons 
are not at the centre of gravity is an additional complication, 
which fortunately can be treated separately as the ' structure 
amplitude ' (James, pp. 49-53, and below, p. 45), and for 
the present it is better to ignore it. In figure 1 the dots 
represent the scattering centres in the plane containing the 
a and b axes of the crystal, and AB represents the direction 
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of a parallel beam of monochromatic X-rays falling on it. 
In general, each scattering centre will diffract rays in all 
directions, the amplitude depending on the relation between 
the direction of the incident ray AB and the direction of 
the diffracted ray BC, but there will be some radiation 
diffracted in every direction. Most of the diffracted radia- 
tion is ' coherent ' with the incident beam, that is the phase 
of the vibration of the diffracted ray BC is fixed by the 
phase of the vibration of the incident ray AB. The total 
amplitude of the radiation diffracted in the direction BC 
will therefore be the sum of the contributions of the in- 
dividual scattering centres, proper attention being paid to 




FlO. 1. — X-rays diffracted by a plane of scattering centres. 

the phase relations. If it happens that for every scattering 
centre B, B', B",— the scattered rays BC, B'C, B"C",-~ 
ar c in phase the radiation diffracted in the direction BC 
*ill be of great intensity, since all the contributions add. 
Un the other hand, if there is no particular relation between 
jhe phases of the rays BC, B'C, B"C",— the radiation will 
^ weak, as the contributions of the various scattering 
centres will mosdy cancel each other out. The condition 
Jjjat the ray BC and the ray B'C should add is clearly 

Bn ^ D + ^ E snou ^ k e an mte g ral multiple of ?,, where 

and BE are perpendicular to A'B' and B'C respec- 

1Vel y. and ;. is the wave-length of the incident X-rays. 
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Our object is to find a simple geometrical way of repre- 
senting the fulfilment of this condition for every scattering 
centre of the crystal. This is best done in stages, a line 
of scattering centres being considered first, then a plane, 
and finally the whole crystal. It is convenient, though 
not necessary, to take the line parallel to one of the crystal 
axes, say a, and the plane including two of the crystal axes, 
say a and b. 

In figure 2, then, OX represents a line of scattering 
centres, at a constant distance a apart. The contributions 
of the scattering centres B and B' will add if the path 
difference B'D + B'E - hX, where h is any integer, positive, 
negative, or zero. (The figure is drawn for h positive ; the 




Fig. 2. — Determination of the path differences between rays 
diffracted by successive scattering centres in a line. 

reader should convince himself that the condition holds for 
other directions of AB and BC. The lines AB, BC and 
OX are not necessarily coplanar.) A little consideration 
shows that if the contributions of the adjacent scattering 
centres B and B' add, the contribution for any other scatter- 
ing centre B" in the row will add also, the path difference 
for B and B" being an integral multiple of that for B and B'. 
The condition that B'D + B'E = hX can be represented as 
follows. From a fixed origin O (figure 3) draw a straight 
line of unit length OP parallel to BA, and from P draw 
a straight line of unit length PQ parallel to BC. Let PP 1 
and QQ' be perpendiculars from P and Q on OX, where 
OX is parallel to the line of scattering centres. Then the 
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triangle OP'P (figure 3) is similar to the triangle B'BD 
(figure 2), since their angles are equal, and so 

OP' _ OP 1, 

Wb " B'B ~ a 

OP' = B'D fa. 
Similarly P'Q' = B'E/a, 

so that OP' + P'Q' - OQ' - {B'D + B'E)/a 

= hX/a 

is the condition that the rays diffracted by all the scattering 
centres B, B', B",'" add. Interpreted geometrically, 




P' Q' 

FlG. 3. — Geometrical representation of the path differences. 

OQ' = hX/a is the condition that Q should lie in a plane 
perpendicular to OX and distant hX/a from O. The 
radiation diffracted by the line of scattering centres has 
therefore its greatest possible intensity if Q lies in one of 
a set of parallel planes perpendicular to the line of scatter- 
ing centres and separated X/a from one another. As Q 
moves away from one of these planes the intensity of the 
diffracted radiation will fall off, and will become inappreci- 
aD je at a distance of the order of X/t from the plane, where 
* ls the length of the line of scattering centres. This is of 
c °urse small compared with X/a, and for many purposes 
ca n be neglected. The falling off in intensity is caused by 
the rays diffracted by the various scattering centres becom- 
mg out of phase with one another, and will be considered 
,n detai l in the next chapter. 
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3. The representation of the condition for reinforcemeni 
that has just been derived is satisfactory for diffraction ol 
X-rays of a single wave-length X, but for some uses th« 
dependence of the spacing of the planes on / is inconvenient 
This may, however, be avoided by a simple device. Tht 
lengths of OP and PO were taken as unity, as any oth 
choice would have seemed artificial as a starting-point, but 
any convenient length may be substituted. If the lengtl 
1 /I is chosen for OP and PQ, OP' becomes B'D/la 
P'Q' becomes B'E/Xa and OQ' = OP' + P'O' becom 
(B'D + B'E)/).a. The condition for reinforcement 
therefore that 00' should be equal to h/a, where h 
an integer. In geometrical language this is the conditio^ 
that O should lie in one of a set of parallel planes per- 
pendicular to the a-axis of the crystal and spaced at distance 
of \/a (figure 4). 




Fig. 4. — Condition that the path differences should be integral 
multiples of the wave-length. 

4. The condition so far derived relates only to a single 
line of scattering centres OX. However, if it is satisfied 
for one line OX, it is satisfied for any parallel line O'X 
(figure 1), since all arguments have involved only th 
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directions of rays, and not the absolute positions of the 
scattering centres. The line OX of figure 1 could be 
translated bodily to the position O'X' without affecting 
the directions of any of the rays in succeeding figures. 
The case is different with regard to lines of scattering centres 
not parallel to OX. Application of the same argument to 
a line of scattering centres parallel to the 6-axis of the 
crystal gives as the condition for the reinforcement of the 
scattered rays that O should lie on one, say the k th , of a 
set of parallel planes perpendicular to the 6-axis and spaced 
a distance i/b apart. If the two conditions for reinforce- 
ment are to be satisfied simultaneously Q must lie on the 
intersection of these two sets of planes, that is on one of 
a set of lines perpendicular to the plane containing a and 
b and spaced as shown in figure 5. If this condition is 







Fiu. 5, — Condition that the path differences between rays dif- 
fracted by all scattering centres in a plane should be integral 
multiples of the wave-length. 

satisfied all centres in the plane shown in figure 1 scatter 
j"- phase, and the diffracted radiation is of maximum 
jntensity. This behaviour is different from that assumed 
*n the derivation of Bragg's law given in chapter I. A 
Kttle consideration makes it clear that reflection as if by 
a mirror corresponds only to the vertical line through the 
or 'gin in figure 5. 
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If all the scattering centres of the crystal are to diffract 
in phase, a similar condition for the reinforcement of the 
rays scattered by a line of centres in the direction of the 
c-axis must also be fulfilled. As now seems natural, the 
condition is that Q should lie on one, say the / lh , of a set 
of parallel planes perpendicular to the f-axis and spaced 
a distance i/c apart. These will intersect the vertical 
lines of figure 5 in a series of points, forming a new lattice 
(figure 6), and all the unit cells of the crystal will scatter 




Fig. 6. — Formation of the reciprocal lattice, expressing the 
condition that the path differences between rays diffracted 
by all scattering centres in a crystal should be integral 
multiples of the wave-length. 



in phase if Q coincides with one of these points. The 
new lattice is characterized by a new set of axes ; that 
having the direction nearest to a is called a*, that nearest 
b is called b*, that nearest c is called c*, and the appropriate 
angles are called **, /?*, y*. The crystal axes and the new 
axes coincide in direction only if the crystal axes are ortho- 
gonal, and their relative magnitudes are generally different. 
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There are, however, a number of interesting and curious 
relations between the lattices, which can be summarized 
in the word reciprocal, and the new lattice, representing 
the conditions that the diffracted intensity should have a 
maximum, is usually called the reciprocal lattice. 

5. The geometrical conditions for strong diffraction are 
summed up in the integers h, k, I, generally called the 
' indices of reflection ', which specify the point of the 
reciprocal lattice with which Q coincides. It is convenient 
to call the direction of the line OQ joining the origin of 
the reciprocal lattice to the point h, k, I the ' hkl direction '. 
The angle between the incident and the diffracted rays 
depends only on the length S of 00, for in the triangle 
OPQ, the angle OPQ is the angle of deviation, and if this 
is called 20, 

S = 00 = OP sin +PQsinO 

= 2 sin 0/X, 
X m (2/5) sin 0. 

If i/5 is replaced by d/n, this is of the same form as 
Bragg's law. The distance 00 should therefore be the 
reciprocal of some interplanar spacing in the crystal 
lattice. 

This, or at any rate the converse of it, that for every 
set of planes with spacing d there are points of the reciprocal 
lattice at distances i/d, 2/d, 2/d, ••• from the origin, is easy 
to show. In constructing the reciprocal lattice we con- 
sidered three non-parallel lines of scattering centres, which 
for simplicity we chose parallel to the crystal axes. If we 
choose them with other directions, the reciprocal lattice 
resulting must be the same, as it expresses the geometrical 
conditions necessary for the reinforcement of the scattered 
rays, and the final edifice must be the same, whatever the 
scaffolding used in constructing it. Suppose then that we 
take the first two lines of scattering centres lying in one 
°f the planes of the set of spacing d, and take the third 
line FG so that one scattering centre belongs tb each plane 
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of the set (figure 7). Let the spacing of the centres in 

this line be d'. The condition for reflection from the 

plane as from a mirror is 
that Q should lie on a 
line ON' through the] 
origin of reciprocal space 
perpendicular to the 
plane, and the condition 
for reinforcement of re- 
flections from successive 
planes is that should 
lie on one of a set of 
parallel planes perpen- 
dicular to the line of 
scattering centres and at 
a distance \/d' apart. 
Thus Q must be one of 
the points Q v Q 2 , Q a , — 
shown in the lower part 
of figure 7, where ON' 
is parallel to FN, OG'\ 
to FG, and the spacing 

of the planes is i/d'. The triangles FJH and OKQi are 

similar, since their angles are equal, so 




Fie. 7. — Construction to show 
that the reciprocal-lattice 
representation includes 
Bragg's law. 



or 



d:d'::i/d':OQ u 
OQ x = i/d. 



The distances of the points Q lt Q 2 , ••• from O are therefore: 
i/d, 2/d, •••, so that there are points of the reciprocal lattice 
at these distances from the origin. These points correspond 
to the first, second, third, ••• orders of reflection from the 
planes of spacing d, but the tendency nowadays is to count 
only first-order reflections, and regard the higher orders 
as being due to different sets of planes, with spacings d/z, 
d/3, '", which happen to be parallel to the set with spacing d. 
Thus for every set of parallel planes in a crystal, there i» 
one point of the reciprocal lattice whose distance from thei 
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origin is the reciprocal of the spacing of the planes, and 
the line joining the point to the origin of the reciprocal 
lattice (the hkl direction as defined above) is perpendicular 
to the planes. The planes are denoted by the same triplet 
of numbers that describes the point of the reciprocal lattice, 
and it can be shown that this triplet is the same as the 
Miller indices of the plane with the order of reflection 
included (James, pp. 5, 51). 

6. The geometrical construction of paragraph 4 deter- 
mines the reciprocal lattice, but for calculations it is 
convenient to have 

the reciprocal axes /t*S 

a*, b*, c* and /' »b* S 

angles a*, /?*, y* 
given explicitly in 
terms of the crystal 
axes a, b, c and 
angles a, £, y. 
The relations be- 
tween the unit 
cells of the lat- 
tices is shown in 
figure 8. Because 
of the method of 
construction of the 
reciprocal lattice 
a* is perpendicu- 
lar to b and c, b* to c and a, and c* to a and b. The 
reciprocal relations follow : a is perpendicular to b* 
and c*, b to c* and a*, and c to a* and b*. Also because 
pf the method of construction, the projection of a* on a 
•s of length i/a, that of b* on b is of length i/b, and of 
£* on c is of length i/c. The reciprocal relations follow ; 
'f the angle between a and a* is ifi, the projection of a* 
°n a is a* cos tjj, so a* cos ifj =• i/a, cos = i/aa*, and the 
projection of a on a* is a cos tp = a/aa* = i/a*. Similarly 
the projection of b on b* is i/b*, and of c on c* is i/c*. 




Fig. 8. — Relations between the axes of 
the crystal lattice and those of the 
reciprocal lattice. Right angles are 
marked. 
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From these relations it is easy to find the lengths of the 
reciprocal axes. Let the area of the face of the parallelepiped 
bounded by b and c be A, that of the face bounded b 
c and a be B, that of the face bounded by a and b be C, 
the volume of the parallelepiped be U, and let the coiv 
responding functions of the reciprocal axes be denoted byi 
the same letters starred. The volume U is equal to the 
area of the face A multiplied by the perpendicular distance 
to the opposite face, and this distance is the projection of 
a on a*, since a* is perpendicular to A. Thus A(i/a*) = U, 

or a* = A/U = be sin tt/U, . . (i) 

and similarly for the other axes. 

If the axes are at right angles each reciprocal axis is 
the reciprocal of the corresponding crystal axis, so that 
U = abc = i/a*b*c* = i/U*. The relation U - i/U* is 
in fact true for non-orthogonal axes, as the following pro- 
cedure shows. Choose rectangular axes with OX along a,, 
OY along b*, and OZ perpendicular to these and with the 
general direction of c and c*. The direction cosines of the 
crystal and reciprocal axes with respect to these axes are 
readily found to be 

a : i, o, o, 

b : cos y, D cosec ft, (cos v. - cos ft cos y) cosec ft, 

c : cos ft, o, sin ft, 

a* : D cosec a, (cos a cos ft - cos y) cosec a cosec ft, 

- D cosec a cot ft, 
b* : o, i, o, 
c* : o, (cos ft cos y - cos a) cosec ft cosec y, 

D cosec ft cosec y, J 
where 
D = \/(i-cos 2 a-cos 2 /3-cos 2 y+2 cosa cos/5 cosy). 

Some of these direction cosines can be written down 
immediately, and the rest follow from known angular 
relations between the axes. The volume of a parallelepiped 
is equal to the product of its edges multiplied by the 
determinant of their direction cosines, so 



(2) 



(3) 
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U = abcD, (4) 

and U* = a*b*c*D 2 cosec a cosec ft cosec y. (5) 

Replacing a* by be sin ol/U, and so on, gives 

U* = a 2 b*c 2 D 2 /U* = i/U, . . (6) 

as it was required to show. Use of this result gives 

a* = A/U = sin ol/oD,-} 

b* =B/U = sin ft/bD,\ . . (7) 

c* = C/U = sin y/cD.j 

From the relation reciprocal to equation 1 

a = A*/U* = b*c* sin a*. abcD 
= a sin a* sin ft sin y/D, 

sin a* = D cosec ft cosec y, . . (8) 

or, since cos 2 a* = i - sin 2 a*, 

cos a* = (cos ft cos y - cos a) cosec ft cosec p, . (9) 

so that 

sin a* = D cosec /5 cosec y, 

sin /i* = D cosec y cosec a, 

sin y* = D cosec a cosec /?, 

cos a* - (cos ft cos y - cos a) cosec ft cosec y, ^(10) 

cos ft* = (cos y cos a - cos ft) cosec y cosec a, 

cos y* = (cos a cos ft - cos y) cosec a cosec /?.- 

7. These results make it possible to write down the 
general expression for the Bragg angle in terms of the 
indices hkl. The length 5 of the line joining the point 
hkl to the origin of the reciprocal lattice may be obtained 
*n the following way. Its components along the reciprocal 
*xes are ha*, kb*, le* ; let it make angles ?,, //., v with these 
^cs. The projection of S on any axis is equal to the sum 
°f the projections of its components, so 

S cos X = ha* + kb* cos y * + le* cos ft*, (11) 

S cos ft = ha* cos y* + kb* + le* cos a*, (12) 

5 cos v = ha* cos ft* + kb* cos a* + le*, . (13) 

an d 5 is equal to the sum of the projections of its 
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components on itself, so 

5 = fia* cos X + kb* cos ft + Ic* cos v. . (14) 

Multiplying (1 1) by ha*, (12) by kb*, (13) by Ic* and adding 
gives 

h 2 a* 2 + k 2 b* 2 + l 2 c* 2 + 2kW*c* cos a* 

+ 2lhc*a* cos /5* + zhka*b* cos y 
= »S(Atf * cos ?. + kb* cos fx + Ic* cos v) 
= S 2 , . . ... . . (15) 

by equation 14. Since S = 2 sin 0/A, the Bragg angle ii 
given by 

4 sin 2 0/A 2 - h 2 a* 2 + £ 2 6* 2 + l 2 c* 2 + zklb*c* cos a* 

+ 2lhc*a* cos/9* + zhka*b* cosy*. . (16 

This equation may be simplified considerably for crystals 
not of the triclinic system. Expressions for each crystal 
system are given in the Internationale Tabellen zur Bestim- 
mung von Kristallstrukturen, and many books on X-ray 
analysis. For the cubic system with a=b>=c,<x.=f}=y 
90 it reduces to 

S 2 = 4 sin 2 O/l 2 = (h 2 + k 2 + I 2 ) /a 2 . . (17) 

8. One reason for the importance of the reciprocal lattice 
is its close relation with X-ray diffraction photographs 
X-ray photographs resemble the reciprocal lattice mucfl 
more than they resemble the crystal lattice ; one coul( 
perhaps say that the photograph is merely a more or lesJ 
distorted reciprocal lattice. For special purposes cameras 
have been designed that reproduce a section of the reciproca 
lattice with little or no distortion (Buerger, 1944). 
general, however, the distortion is more or less severeM 
and before proceeding with the development of thu 
reciprocal lattice it may be well to describe briefly certain 
types of X-ray camera and the photographs they produce 

9. In the simplest possible arrangement the crystal an' 
the film are both stationary and a narrow beam of practically 
parallel X-rays falls on the crystal. If the beam contains 
a wide range of wave-lengths this gives a so-called Lau 
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photograph (James, p. 19), which is not nowadays of 
primary importance. If the beam is monochromatic and 
the crystal perfect there will be in general n<> strong diffrac- 
tion maxima. The reciprocal lattice and the line OP 
(paragraph 3) are fixed, and strong diffraction would occur 
only if some point of the reciprocal lattice were at a distance 
i/A from P, or in other words were to lie on a sphere 
(called the ' sphere of reflection ') with centre at P and 
passing through the origin of the reciprocal lattice. For 
an arbitrary orientation of the crystal this is an unlikely 
coincidence. If the crystal is not perfect, however, the 
region of appreciable scattering is not concentrated in a 
point of the reciprocal lattice, but extends to a greater or 
less distance from it, and it is quite likely that the sphere 
of reflection will intersect one or more scattering regions, 
and give rise to more or less diffuse streaks or spots on 
the photograph. Cameras specially designed for investi- 
gating such diffuse regions in reciprocal space are now in 
use, embodying ingenious devices for obtaining a beam as 
nearly monochromatic as possible. Any continuous radia- 
tion would produce Laue spots and make the interpretation 
more difficult. 

10. If the crystal is oscillated or rotated while the ex- 
posure is being made the reciprocal lattice will oscillate or 
rotate through the same angle, many points of the reciprocal 
lattice will pass in turn through the sphere of reflection, 
and the strong scattering occurring as each point intersects 
the sphere will produce corresponding points on the film. 
In the commonest type of moving-crystal camera the crystal 
ls mounted in the centre of a cylindrical * film, adjusted 
s o that one axis is vertical, and rotated about this axis 
during exposure to a horizontal beam of X-rays. An 
0s cillation with uniform angular velocity over a convenient 
f ange (5 , i ° or 15 ) may be substituted for a complete 
Station. Suppose for definiteness that the c axis is vertical . 

A flat film may be used, but the ' layer-lines ' are then curved 
Uames, pp . 2 ^y 
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Then the plane of the reciprocal lattice containing a* ancfl 
b* (the plane / = o) is horizontal, and OP lies in it (figure 9) J 

The sphere of reflection! 
cuts it in a great circle 
and cuts the parallel plana] 
for which / = ± 1, ±2, 
in circles of progressively 
diminishing size. Thi 
diffracted rays for / = o all 
lie in a horizontal plane. 
and strike the film alongB 

the 'equator'. The 
— ^- diffracted rays foij 
other values of / lie 
on right circular cones 
whose angle depends on 
the values of / and c ; tin 
semi-angle is easily seen tg| 
be arc cos (ll/c). Thea 
cones intersect the film ii 
circles whose planes artfl 
parallel to the equatori 
plane, so that when thfl 
film is opened out th 
diffraction spots lie on 1 
series of parallel ' laye 
lines '. The length of th 
c axis is readily obtained 
from the spacing of thes 
lines and the radius of th 
camera. The / index of 
spot is equal to the nu 
ber, counted from tfafl 
equator as zero, of the line 




Fig. 9. — Formation of layer lines 
on a moving-crystal photo- 
graph. The X-rays emerge 
from the sphere of reflection 
at O, P is the centre of the 
sphere, and four possible 
positions of Q are shown, 
leading to spots with 
1 = 2, 1, o and — 2. As the 
crystal moves the reciprocal 
lattice rotates with it, bring- 
ing other points to the surface 
of the sphere of reflection, 
and producing other spots on 
the layer lines. 

on which it lies. Thfl 

other indices of the spots are generally obtained by 

method described by Bernal (1926) which involves 
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certain amount of trial and error, but a direct method 
has been devised (Bunn, Peiser and Turner- Jones 1944 ; 
Turner-Jones and Bunn 1946). 

11. In the Weissenberg camera a cylindrical screen is 
placed between the crystal and the film, with a slot placed 
so as to allow only the spots of one layer line to pass 
through. The film is moved parallel to its axis as the 
crystal is rotated, so that the spots are all over the film, 
instead of being concentrated on a line. It can be shown 
that they lie on two families of curves, from which the 
pair of indices not determined by the layer line can be 
found directly. The method is used chiefly in structure 
determinations, and will not be further described here 
(Weissenberg, 1924). 

12. The methods mentioned so far require a fairly large 
crystal (say not less than o-oi cm. in linear dimensions), 
and if this is not available the powder method can be used. 
In this the specimen, consisting of a large number of small 
crystals of random orientation, is placed in the centre of 
a cylindrical film and bathed in a narrow beam of approxi- 
mately parallel X-rays. During exposure the specimen is 
rotated to make the orientation still more random. By 
this means nearly all directional effects are averaged out, 
the angle of deviation being the only remaining variable. 
The intensity of the X-rays falling on the film with a 
deviation between 20 and 2(0 + dO) will depend on all the 
"ttensity that, in the representation in reciprocal space, 
lies between spheres of radii 5 = 2 sin 0/X and 

S + dS = 2 sin e/k + d(2 sin B/X). 

" 1(0) dd is the intensity diffracted with a Bragg angle 
between and 6 + dO, and I(S) dS is the intensity between 
concentric spheres with radii S and S + dS, then except 
to r certain geometrical factors, 

7(0) dO = 7(5) dS 

= (2 cos d/?.)I(S)dd, 
So that 1(0) =(2 cos 0/k)I(S). . . (18) 
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In favourable cases the spheres will not intersect more thanj 
two regions of high intensity in reciprocal space (say thosej 
near h, k, I and near - h, - k, - 1), and the photographs 
will consist of a series of lines. Usually the regions round 
h, k, I and round -h, - k, -I are related by a centre of 
symmetry, and their superposition on the film causes no 
trouble. In unfavourable cases the spheres will intersect 
regions corresponding to several values of h, k, I and the 
interpretation of the film becomes more complicated, but 
powder photographs are of some use as long as discrete 
lines are visible, even when the* lines are due to the super- 
position of several regions of high intensity. 

For purposes of structure analysis the important quantity 
connected with a line is its total intensity, that is I=\I(0) dO 
the integral being taken over a sufficient length of film to 
include the whole line. Except for certain geometrical 
factors, this is the intensity that would be scattered by 
one unit cell at the angle multiplied by the number of 
unit cells in the specimen. We shall be more concerned 
with the shape of the line, 1(0), and especially with its 
breadth, or the range of angle over which the reflection 
extends. Various definitions of this have been used, the 
earliest being the half breadth. If I is the maximum 
value of 1(0), and V 2 are the values of for which 
1(0) = U , the half breadth is defined as 

/ 3'=2(O a -0 1 )- • • • (»« 
A definition of breadth better adapted to calculation is the 
integral breadth, defined by 

/*=2/// . . . • (20) 

(The factor 2 occurs in these expressions because they were 
originally defined in terms of the angle of deviation, not 
the Bragg angle.) If I„ is the maximum value of I(S), 

/? = /A//, cos 0, . . . (21) 
so that the apparent particle size is 

e = X/fi cos = /,//. . . . (22) 
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13. In deriving the relation between 1(0) and I(S) above 
it was tacitly assumed that the variation of intensity on the 
film depended only on the distribution of intensity in 
reciprocal space. This is never true in practice ; the lack 
of parallelism in the incident beam, the finite length and 
diameter of the specimen, and the range of wave-length 
included in even a crystal -reflected X-ray beam would lead 
to a line of finite breadth even if the specimen were ideal. 
The line profile obtained by photometering the film differs 
from that which the crystal imperfection alone would 
produce, and it is necessary to correct it before using it to 
obtain information about the specimen. Several methods 
have been devised for doing this, depending on the use of 
two lines, one produced by the specimen under investiga- 
tion, and the other produced under the same experimental 
conditions by a specimen as nearly ideal as possible. From 
the two line profiles it is possible to derive the line profile 
due to specimen imperfection alone, but the process is 
tedious and one has generally to be content with the integral 
breadth only. It is found that the breadth of a line pro- 
duced by a specimen with imperfections is not simply the 
sum of the breadth due to experimental conditions and the 
breadth due to specimen imperfection. The squares of the 
breadths are more nearly additive, but the exact relation 
depends on the type of imperfection and the experimental 
arrangement. For routine work the correction procedure 
devised by Jones (1938) is perhaps the best. 
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CHAPTER III 

RECIPROCAL SPACE FOR SMALL CRYSTALS 

i. As was seen in the preceding chapter, the radiation 
diffracted by a crystal is at a maximum when the relation 
between the directions of the incident and diffracted rayj 
is such that the point Q, found by drawing a line OP oi 
length x/7. from the origin of the reciprocal lattice parallel 
to the incident ray reversed, and a line PO from its terminus 
parallel to the diffracted ray, coincides with one of the 
points of the reciprocal lattice. If the coincidence is not 
exact there is still some diffracted radiation, and the next 
step in the investigation is to find the way in which the 
intensity of the diffracted radiation depends on the position 
of O. Suppose that Q is nearly at the point h, k, I of the 
reciprocal lattice, its co-ordinates in terms of the reciprocal 
axes being h + u, k + v, I + w. The successive scattering 
centres in the line OX of figure 2, chapter II, will no longer 
be scattering exactly in phase, and the bigger // is, that 
is the further Q departs from the plane for which 
B'D + B'E = hX, the greater will be the phase difference 
between the rays scattered by successive centres. I 
u were equal to unity, Q would lie on the plane correspond- 
ing to h + 1, and the phase difference would be one wave- 
length, or 27T radians. For smaller values of u the phase 
difference between rays scattered by successive centres 
will be proportionately less, and so equal to 27111 radians- 
If the length of the line of centres is T the number of them 
is T/a, and the amplitude of the scattered radiation is thl 
resultant of N = T/a equal vibrations each 6 ■ zna radian* 
out of phase with the preceding one. 

2. A graphical method of finding the resultant G of M 
equal vibrations with successive phase differences 6 is illuS* 
trated in figure 1. Each vibration is represented by a linl 
of unit length, and the phase differences are represented ty 
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Fig. 1. — Summation of JV 
out-of-phase vibrations. 



RECIPROCAL SPACE FOR SMALL CRYSTALS 25 

drawing each line at an angle d to the preceding one. The 
resultant amplitude is given by the distance between the 
free ends of the first and 
last line. In the present 
case N is large, and the 
polygon made from the 
individual vibrations may- 
be considered as the arc of 
a circle, at any rate as long 
as 11 is small compared with 
unity. The angle subtended 
by the arc at its centre is 
Nd, and its length is N, so 
the radius of the circle 
is i/d. The length of the secant joining its ends, which 
is the amplitude of the resultant, is twice the radius 
multiplied by the sine of half the angle subtended at the 
centre, or 

2 . Nd. 

r m_ . . . . (I ) 

More accurately, the angle subtended by the arc of the 
circle is Nd, but its length is AT.2Rarcsin (1/2R), where 
R is the radius of the circle. The radius is therefore 
given by 

zNR arc sin (1/2R) = N3R, 
0r R = 1/2 sin (6/2). 

fhe exact expression for the reflected amplitude is thus 
sm(Nd/z) 
sin (d/2) 

b ut the expression (i) is a sufficient approximation for all 
a Pplications made here. 

3- On substituting the values given above for N and d 
the resultant amplitude becomes 

sin (7iTu/a) t 
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(2) 
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The reflected intensity is proportional to the square o! 
this, that is to sin a faTu/a) 



{nuy 



(3) 



For u = o this has the value {T/a) 2 , and as u increases it 
rapidly decreases, reaching zero for u = a/T. It then | 
increases a little, but never amounts to more than five pel 
cent of its value for u = o. For large u it falls to zero 
IT"*. Figure 2 shows its behaviour as a function of //. 




Fig. 2. — Resultant of N = T/a out-of-phase vibrations as 
a function of u. 

4. We have so far considered the behaviour of a single 
line of scattering centres. In order to investigate thfi 
behaviour of a complete crystal it is necessary to know i 
shape. In principle it is/ possible to work out the scatte 
ing from a crystal of any shape, and it has indeed been do 
for many (Patterson, 1939). To illustrate the proble 
two shapes will be considered here, the parallelepiped and 
the sphere. Others can be worked out on the same prin- 
ciples, but the mathematics becomes complicated. The 
general expression is derived in section 6 below. 

Suppose, then, that the crystal has the form of a par- 
allelepiped, whose edges are parallel to the crystal axes, 
and are of lengths X, Y, Z. The value of T for all lines 
of scattering centres parallel to the a-axis will be X, an< 
the dependence of the amplitude of reflection on u will b 

sin (nXu/a) 
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for each such line. The dependence of the amplitude on 
v is found by adding Y/b vibrations, each of amplitude 
sin (nXu/a)/7iu, with a phase difference of 2tzv between 
successive vibrations. By the construction of section 2 
the resultant is 

sin (nXu/a) sin (nYv/b) 

Till TtV ' ^' 

Finally, the dependence of the amplitude on w is found by 
adding Z/c vibrations, each of amplitude given by expres- 
sion (5), with a phase difference of 27izv between successive 
vibrations. The dependence of the diffracted amplitude 
G on u, v, w is therefore 

G(u, v, w) = Sln fog*/*) sin {nYv/ b) sin (nZw/c) 
nu nv nw ' 

The diffracted intensity as a function of position in 
reciprocal space, H(u, v, w), is proportional to the square 
of this, that is 

H(«, v, w) = Sin2 (**»/«) sin2 faYv/b) sin 2 (nZ w/c) 
{nu)* (nv) 2 (7m) 2 
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" 



(7) 

The maximum value of H(u, v, w), for u = v = w = o, is 
\XYZ) 2 /(abc) 2 . Since XYZ is proportional to the volume 
v of the crystal, and abc to the volume U of the unit cell, 
the maximum value of H is the square of the total number 
°f unit cells in the crystal. For a given distance from the 
Point of the reciprocal lattice H has its greatest values 
along the reciprocal axes, that is, the intensity falls off least 
r aptdly in directions perpendicular to the faces of the 
crystal. This is found to be true in general ; each plane 
. a ce of a crystal is represented by a ' spike ' of increased 
jntensity in reciprocal space, the direction of the spike 
e 'ng perpendicular to the face (von Laue, 1936 ; Ewald, 
'938, 1940). If the crystal is plate-like, with, say, X and 
^ much greater than Z, the intensity of reflection as 
^ unction of position in reciprocal space is practically 
°ncentrated on the c*-axis, forming a line perpendicular 
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to the plate and of length inversely proportional to its 
thickness. If the crystal is needle-like, with, say, X and Y 
very much less than Z, the intensity of reflection is con-' 
centrated into a plate perpendicular to c*, extending a 
little bit further in the directions of a* and b* than in the 

intermediate directions. 
These plate-needle and 
needle-plate relations are J 
a further illustration of 
the reciprocal nature of] 
the relations between the] 
two sets of axes. 

5. For a spherica 
crystal the calculation i 
a little more complex,] 
and to simplify it as far! 
as possible it will be| 
assumed that the axes of the crystal are orthogonal. Sup 




Fig. 3. — Calculation of the in- 
tensity of reflection from a 
spherical crystal. 






pose that the diameter of the sphere is D. The length o 
a line of scattering centres at a distance r = V(Z)/2) 2 - *| 
from the a-axis is 2x (figure 3), so the amplitude diffracte 
by this line is, as a function of u, 
sin {zTixu/a) 



(8) 



Till 



The amplitude diffracted by the whole crystal, fo| 
v = iv = o, is got by summing this expression for all lines 
of scattering centres. The number of lines of scattering 
centres passing through unit area perpendicular to a is 
equal to a/U, so that the sum is approximately 



-r 

J D/2 



<D/ 

= {zaf 



sin (zTixu/a) a . 

— J — ~27trdr 

71U U 



'uU) f° sin {2nxu/a)V'(D/2¥ - * 3 dV(D/z) 2 ^*1 

J Dl'l 



( 2 a/«t/)j 



■D/ 
•X>/2 



x sin (27txu/a) dx. 



(9) 
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The integration is easily performed by parts, and gives as 
the amplitude of the diffracted radiation 

n . . 71 sin (nDu/a) - (tcDu/o) cos (nDu/a) . . 
G{U) =W- ~Waf * (I0) 

This has been derived as the amplitude diffracted as a 
function of u for v = w = o. For a spherical crystal, 
however, there is nothing special about the a-axis, and the 
amplitude of the diffracted radiation will depend only on 
the distance s = u/a of O from the point of the reciprocal 
lattice. This amplitude G(s) will therefore be given by 
equation (10) with u/a replaced by s. The maximum 
value of G, being the limit of equation (10) for u = o, is 
7iD s /6U, which is, as would be expected, just the total 
number of unit cells in the crystal. At the points of the 
reciprocal lattice all unit cells scatter in phase, whatever 
the crystal shape, and the intensity of reflection is propor- 
tional to the square of the number of unit cells. The total 
intensity of a reflection (that is, the integral of the intensity 
over all the region of reciprocal space for which the 
intensity is appreciable) is proportional to the actual 
number of unit cells, so the ' integral breadth ' in reciprocal 
space of a reflection is N/N' 2 = i/7V, whatever the crystal 
shape. The shapes of the regions of high intensity, and 
hence the details of the X-ray photographs, are of course 
dependent on the crystal shape. 

6. Though it is not needed for applications in future 
chapters it is perhaps worth while to write down the general 
expression for the amplitude of reflection as a function of 
Position in reciprocal space. For simplicity the treat- 
ment will be limited to orthogonal axes. The rays scat- 
tered by the unit cell that is »j unit cells in the a direction, 
n 2 in the b direction, and w 3 in the c direction away from the 
Ce 'l at the centre of gravity of the crystal will be advanced 
ln phase with respect to those scattered by the cell at the 
Ce ntre of gravity by 



<5„ = 271(71^1 + n 2 v + n 3 w) 



(11) 
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radians. The amplitude scattered by the whole crystal 
will be the sum of N = V/U contributions (where V is the 
volume of the crystal and U the volume of a unit cell) of 
equal magnitudes, and phases given by equation (n) 
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Fig. 4. — Position of the n" 1 unit cell. 



These cannot be added by the construction given in para- 
graph 2, since d„ depends on three variables, n x , n 2 , n 9f 
but the addition can be done trigonometrically by noting 
that the resultant of a number of lines of lengths F v 
F 2 , ••• at angles of d v d 2 , ••• with the horizontal is repre- 
sented in magnitude G by 

C n -=A*- + B°- . (1 J 

and in phase 6 by 

tan <5 = B/A (1 J 

where A = F, cosd x + F 2 cos<3 a + ••• . (14) 

and B = F 1 sin d x + F 2 sin <5 2 + •••. 

In the present case the lengths are all unity, and the phases 
are given by equation (11), so 

A = the sum of cos 27i(n x u + n 2 v + n z w) for all values 

of n v n 2 , >/ 3 , 
and B = the sum of sin znfau + n 2 v + n 3 w) for all values 

of n x , n 2 , n 3 . 
If the crystal is reasonably large u, v, w are small compared 
with unity, and the summation may be replaced by integra- 
tion without appreciable error. The number of unit cell* 
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in the volume element dx dy dz is dx dy dz/U, and x = n x a, 
V - n 2 b, z = n 3 c, so 

A= U */// cos 27i(ux/a+vy/b+wz/c) dx dy dz, (15) 
B = U- 1 }]] sin 27i{ux/a +vy/b +wz/c) dx dy dz, ( 1 6) 

H(u,v,w)=A*+B*, ( I7 ) 

where the integrations extend throughout the volume of 
the crystal. If the shape of the crystal has a centre of 




Fig. 5. — Summation of the contributions from N unit cells. 



symmetry B vanishes. It will be noticed that if u, v, w 
are replaced by - u, - v, - w, A is unaffected and B is 
merely changed in sign. II(u, v, w) has therefore a centre 
of symmetry at the point of the reciprocal lattice, even if 
the crystal shape is not centrosymmetric. The mathe- 
matically-inclined reader may verify that equations (15)- 
(17) lead to (8) and (10) for crystals having the shapes of 
a parallelepiped and sphere respectively. 

7. It will be noticed that none of the expressions giving 
the amplitude or intensity as a function of u, v, w depends 
explicitly on h, k, I. The surroundings of each point of 
the reciprocal lattice are therefore identical, except perhaps 
tor a factor depending on the ' structure amplitude \ This 
should be noted, as it contrasts with the appearance of the 
reciprocal lattice when the crystal is distorted or contains 
mistakes '. 

REFERENCES 

• P. Ewald, 1938. Annates tie I'Instittit Henri Poincart, VIII, 79. 
7~ 1940- Proc. Phys. Soc, 52, 167. 
2 1 - von Laue, 1936. Ann. Phys. Lpz., 26, 55. 
A - L. Patterson, 1939. Phys. Rev., 56, 972. 



POWDER PHOTOGRAPHS OF SMALL CRYSTALS 



33 



CHAPTER IV 

POWDER PHOTOGRAPHS OF SMALL CRYSTALS! 

i. It has been seen that X-ray diffraction is concentrated 
into small regions of reciprocal space in the immediate 
neighbourhood of the points of the reciprocal lattice. The 
diameter of such a region is of the order of i/r, where t is 
a linear dimension of the crystal, and the distance between 
the points of the lattice is of the order of i/a, where a is 
a lattice parameter. The crystals normally used in moving- 
crystal cameras are not smaller than io~ 2 cms., and a is 
not normally larger than 2 x io -7 cms., so the regions of 
strong diffraction have a magnitude of less than io -4 of 
their distance apart. The region of strong diffraction is 
therefore inappreciably small, and the finite size of the 
spots on rotating-crystal photographs is not due to diffrac- 
tion in the ordinary use of the word. Geometrical optics 
suffices to explain the size and shape of the spots. With 
powder photographs, however, there is no lower limit to 
the size of crystals that may be used, and if t is of the 
order of io -5 cms. the size of the regions of large intensity 
becomes appreciable in comparison with their distance] 
apart, and the lines on the powder pattern are broadened.1 
It is necessary therefore to determine how the broadening 
depends on the size and shape of the crystals producing 
the powder pattern. 

2. The intensity diffracted with a Bragg angle between 
and 6 + dO on a powder photograph is proportional to 
the intensity that, in the representation in reciprocal space, 
lies between concentric spheres of radii 5 = 2 sin 6/?. and 
S + dS = 2 sin 6/). + (2 cos 0/X) dO, and the relation be- 
tween the line-profile 1(0) dd and the total intensity included 
between the spheres I(S) dS is (chapter II, paragraph 12) | 

1(8) = (2 cos 6 /X)I(S). . . (1) 
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The line-profile on the powder photograph can be got 
therefore if the integral of H(u, v, w) over the space between 
the spheres can be determined. With expressions for H 
as complicated as equation (7) or the square of equation (10) 
of chapter III this is difficult, but it has been done for 
the second one (Patterson, 1939). Normally, however, 
S is very large compared with the regions over which 
// is appreciable, and no error is introduced by replacing 
the spheres K lf K 2 by their tangent planes L v L 2 (figure 1). 




Fig. i. — Determination of the intensity of reflection as a function 
of from the intensity as a function of position in reciprocal 
space. 

The integral over the region between the planes is difficult 
enough, and is in fact practicable only if the crystals are 
spherical, or for certain reflections from a parallelepiped. 
For spherical crystals it has been done by Patterson (1939), 
and attempts for other crystal shapes have been made by 
replacing the correct expressions for H by approximate 
ones that arc easier to integrate (von Laue, 1936 ; Jones, 
1938 ; etc.). This procedure is, however, of doubtful 
validity, and gives results that are correct only in order of 
magnitude. 

The direct calculation for spherical crystals is somewhat 
too complex to reproduce here, but it gives the same result 
as we shall obtain by a simpler method. Instead of cal- 
culating the amplitude of reflection from the whole crystal 
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as a function of position in reciprocal space, squaring to 
get the intensity, and integrating to get the intensity with 
S between S and 5 + dS, we shall take the intensity that 
is scattered with S between S and S + dS by a single 
line of scattering centres in the crystal, and sum it for all 
parallel lines. It is obvious that the second process con- 
tains fewer mathematical steps, and it happens in addition 
that they are easier to carry out. 

3. The problem is to calculate the integral breadth of 
a general reflection (that is, one with indices h, k, I) from 
a crystal of arbitrary shape. The amplitude diffracted by* 
a line of scattering centres has been found to be 



G(u) 



sin {nTu/a) 



mi 



W 



where T is the length of the line, a is the repeat distance 
in the direction of the line, and u is the fraction of the 
distance i/a between parallel planes perpendicular to the 
line by which O misses the nearest point of the reciprocal 
lattice. The intensity scattered with u between u and 
u + du is then 



//(;/) du 



sin 2 (nTu/a) 
(?t«) 2 



du. 



(3) 



Suppose that the line of scattering centres is parallel to 
the line joining the origin of the reciprocal lattice to the 
point h, k, I, and that x, y, z are orthogonal axes in the 
crystal with x in the direction of the line. The total 
intensity scattered by the crystal with u between u and 
u + du will be got by summing expressions like equation (3) 
for all such lines of scattering centres in the crystal, or to 
a sufficient approximation by multiplying (3) by the 
number of lines per unit area and integrating with respect 
to y and z throughout the crystal. The number of lines 
of scattering centres passing through unit area perpendicular 
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to the lines is a/U, where U is the volume of the unit 
cell, so 

ri \ j a f f s ' n2 OrtTu/a) , , , 

m *"v\) W * i ** h ■ (+) 

This quantity I(u) is very closely related to I(S) ; in fact, 
the only difference is that u is expressed as a fraction of 
the distance i/a between parallel planes, and S is an actual 
distance in reciprocal space. By replacement of du by 
a dS and u/a by s equation (4) becomes 

m 4SrU-*lpgp**i9, . (s) 

s being the actual distance by which S misses the point 
of the reciprocal lattice.* 

The maximum intensity I s is given by equation (5) with 
s = o, so 



I 8 'U-^T'dydz, 



or, since jdx = T, 



I s = U-'tf$Tdxdydz 



u- 



TdV, 



(6) 



(7) 



where the integration extends throughout the volume of 
the crystal. The total intensity / of a reflection is equal 
to the number of unit cells in the crystal, V/U, so that the 
apparent particle size is 



e=I g /I= V-H TdV. 



v 



(8) 



* In this paragraph it has been tacitly assumed that it is possible 
l o pick out a line of scattering centres having the direction of the 
line joining the origin of the reciprocal lattice to the point h, k, I. 
fror a cubic crystal this can be done accurately, and for other 
types of crystal to any desired degree of approximation. If 
Co ,uation (4) holds for directions differing slightly from that 
desired, but enclosing it on all sides, considerations of continuity 
requi rc that it hold for the correct direction also. For a derivation 
ln which this difficulty does not arise, see Stokes and Wilson 
('944). 
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The physical significance of (8) is readily seen : the* 
apparent particle size is the volume average of the thick-] 
ness of the crystal measured perpendicular to the reflecting) 
planes. This is a very elegant result, and carries con- 
viction with it, since it seems intuitively probable that the 
broadening of X-ray reflections due to the smallness of the 
diffracting crystal will depend only on the form andj 
orientation of the crystal, and not on its particular internal 
construction. The form 

6 = V-^jT 2 dy dz . . . ( 9 )| 

obtained from (6) instead of (7), or that given in equation | 
(14) below, is frequently more convenient to evaluate. 
The maximum intensity as a function of 6 is 



1'0\ 

e, oft 



(10) 



(") 



I = (2 cos B/X)I S 

= (2cos0/AC/)f TdV, 

and the integral breadth of the reflection is 

P = 2///0 
IV 

cos of TdV 

4. The apparent particle size is easily calculated fo: 
simple geometrical shapes. For example, for a sphere it 
will be the same for all reflections, and given by (figure 2) 

6 = J^/JT 2 dy dz = (6/jiZ) 3 ) | [ W \d o - - ^)r dr d6 

Jo J 

J2T.D 4 _D 4 " 

D*18 76 

= 3-0/4, • . (ial 

where D is the diameter of the sphere. It is usual, how-} 
ever, to take p m tyV as the standard ' linear dimension '] 
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the crystal, and in terms of this 



(13) 



= -3 -P 
4\7 71 



= 0-930/). 

The apparent particle size e is thus just a little less than 
the actual particle size p. 

For crystals other than spherical e will depend on the 
indices of reflection h, k, I. For reflection from a cube 
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FlG. 2. — Calculation of the apparent particle size for reflections 
from spherical crystals and for the no reflections from 
cubic crystals. 

face it is clearly p, since p is the cube edge, and for planes 
perpendicular to a face diagonal it is (figure 2) 



e = 2p- 



I) 



(\/2p - 2y) 2 dy dz 



(H) 



= Wzp/3 

= o-943/' 

For this reflection also the apparent paiticle size is some- 
what less than the actual. 

5. For the practical calculation of the factor connecting 
£ and p for a general hkl reflection the form of equation (9) 
is not so convenient as one that can be obtained from it 
by a simple transformation. In figure 3, let A be the 
Cr ystal, and let B be its ' ghost ', shifted a distance t in 
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Fig. 3. — Crystal and its ' ghost 

the hkl direction. Let V{t) be the volume, shaded in the 
figure, common to the crystal and its ghost. Then 

V{t)=lST'dydz 

= SKT~t)dydz, . . (J 

where T is the length of a line drawn through the point 
x, y, z in the hkl direction and bounded by the crystal 
and the ghost. For t = o, V(t) is the volume of the 
crystal, and V(t) vanishes for a value of t depending on 
the size and shape of the crystal. Suppose now that 
equation (15) is integrated with respect to t from t = o 
to a sufficiently large value to make V{t) vanish. This 
may be symbolized by 

+ V{t)dt. 




dydz. . (16) 



f 



Then [V(r) dt = j IT7T ' dt - p dt 

Clearly dt is simply T, the thickness of the crystal 
Jo 

measured in the hkl direction through the point x, y, z f \ 
and j 'tdtis T z /2, so 

f V{t)dt = \[['r-dydz 

= Ve/2. 
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If the ghost is shifted a distance t in the opposite direction 
the common volume is the same as before, so the integral 
of V(t) from a sufficiently large negative value to t = o is 
also Ve/2. Thus 



[°V(l)dt+i + V(t)dt, 



or for shortness $V(t) dt, is equal to Ve, and 
6 = v~ l $V{t) dt. . 



■ (18) 



This expression may be used to evaluate the relation 
between s and p for the ///;/ reflection from a cube. Sup- 
pose that the indices are positive and arranged in the order 
h >k > I. This is merely a matter of naming the axes. 
The volume common to the crystal and its ' ghost ' is 
then a rectangular parallelepiped of edges p - h' \t\, 
I' - k' \t\, p - I' \t\, where /»', k', V are the direction cosines 
of the reflection, being //, k, I divided by -\/(h 2 + A 2 + / 2 ). 
The limits of integration are t = ± p/h', and the apparent 
particle size is therefore 

J+p/h- 
(p-h'\t\)(p-k'\t\)(p-l'\t\)dt 
-p/h- 

= 2 [pi/h'-ih'+k'+l'W^h't+ih'k'+k'l'+l'h'W/ih'* 
/> 3 L 

-h'k'iy/4h'*] 

__ V(h* + k°- +¥){6h i - 2hk + kl - zlh) 



6// 3 



>. • • (i9) 



Table I gives values of the ratio p/e for a number of 
reflections from crystals whose external form is cubic, 
octahedral or tetrahedral. It is somewhat modified and 
extended from that given by Stokes and Wilson (1942), 
^vhere explicit formulae for these shapes and also for the 
Octangular parallelepiped are given. 
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TABLE I 

The ratio of the true to the apparent particle size for various reflection! 
from crystals having the form of a tetraliedron, a cube or 
octahedron. For a sphere the ratio is constant at i -0747 

Reflection Tetrahedron Cube Octahedron 

100 1*3867 10000 1-1006 

no 0-9806 1-0607 1-0376 

hi 1-2009 1-1547 1-1438 

210 1-2403 1-0733 i'i°75 

211 11323 1-1527 I-IOOI 
221 1-1556 1-1429 1-1185 
310 13156 10672 11138 
3U 12543 11359 11211 

320 1-1538 1-0698 10902 

321 1-1119 11394 i'°955 

410 1-3453 1-0583 1-1123 

322 1-1772 1-155° 1-1304 

411 1 3074 i-i 174 1-1207 

331 I- ii35 1-1262 10963 
421 1-2104 1-1324 I-II33 

332 11826 I'I5I3 1-1334 
430 1-1094 1-0667 1-0786 
43i 10878 1-1240 10835 
510 1-3597 10506 i-iioi 

6. The line-profile, that is the intensity of reflection 
a function of the angle of reflection, is given by equation I 

1(6) = ( 2 COS0/x)/(S), 

where I(S) is the integral of the intensity of reflection 
reciprocal space over the region bounded by spheres of 
radii S and 6' + dS. Equation 5 gives a general expression 
for I{S), 



(tzTs) 



jffsin 2 (:-r 

U)) (Ttsf 



dy dz, 



but, like the corresponding equation (9) for the apparent 
particle size, this expression is not usually the most con- 
venient to evaluate. One alternative form, corresponding 
to equation (7), is 



I(S) 



u 



rpS^&Ts) 

v (nTs)* 



(20) 
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where the integration extends throughout the volume of 
the crystal, but the most useful form is probably 

I{S) = U- l jV{t) COS (27ZSt) dt, . . (21) 

where V(t) is the volume common to the crystal and its 
' ghost ' shifted a distance t in the hkl direction, and the 
limits of integration are the positive and negative values 
of / for which V(t) vanishes. This may be shown by the 
same method as was used in the preceding paragraph. 
Since 

V{t)-U(T-t)dyd* 

(equation 15), 

j V(t) cos (271st) dt = j J J (T - t) cos (271st) dt dy dz 



T 1 cos (271st) dt - I t cos (271st) dt dy dz 
J 11 Jo 

= jT f ^m^sT) _ ^mi^tsT) + r+ sin (271st) ^\ ^ 

JJL 271S 271S Jo 271S J 

-'ff L 



- cos (ztzsT) 



dy dz 



(22) 



dydz 
(my 



^rrsin^.r) 
J J (*0 2 y 

Since I V(t) cos (271st) dt is equal to I V(t) cos (znst) dt, 

i (F(0cos(2,^)^ = JJ sin2(?Kr) 

= UI(S), . . (23) 

as it was required to show. 
< 7- Little attention has so far been given to line-profiles, 
s mce integral breadths are easier to investigate experi- 
mentally. Equation (21) may, however, be evaluated for 
any given value of V(t). Let t be the value of t for which 
4 
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V(t) vanishes. Then, by repeated integration by parts, 
I(S) = U- x jV(t) cos (wist) dt 

= 2t/- 1 fV(*)C0S(2JW*)<fc 

Jo 

= 2t/-l[W in {2mt) + V'{tf™^ 
L 271S W (27CS) 2 



_ K » (f) sin (271st) _ 



(2.T5) : 



'(/) 



cos (271st) 
(27a)* 



r v . lv{t) cos(2^) dt 



where the primes denote differentiation with respect to t.\ 
The first term in the square brackets vanishes at both 
limits, since V(r) is zero. For ordinary crystal shapes, 
such as all those mentioned in paragraph 5, V(t) is a cubic 
in t, so that V" is a constant, V lv (t) is identically zero 
and the final integral vanishes. The line-profile is therefore 
given by 

I(S) = \U-*[- (ns)-W(o) + (tk)-2 cos (27isr)V'(r) 
-K^)- 3 sin (27isr)V"(T)+\(ns)-*{i -cos (zjin)}V'"]. (25) 

It is easy to find I(S) for any particular case by substituting 
the appropriate values for V(t) and x in equation (25) 
For example, if the crystal is a cube of edge p, 

V(t)=(p-h't)(p-k't)(p-l't) 
and r = p/h' (paragraph 5), so 
V'(t) = - h'(p - k't)(p - Vt) - k'(p - l't)(p - h't) 

- l'(p - h't)(p - k't 

V'(o) = - p2(h' + k' + l'), 

V'(r) = - p*(h' - k')(h' -l')/h\ 

V"(t) = 2h'k'(p - Vt) + 2k'l'(p - h't) + 2l'h'(p - k't), 

V"(r) = zp(h'k' - 2k'V + l'h'), 
V" = - bh'kT. 
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The line-profile is therefore 
P 2 
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W 



±VJ1_ 

2U\_(71S)' 



-QcOs(27ips/h') 



(7U) 



!,' 



■Ui'k' - 2k'V + l'h') sin (znps/h') 



jh'k'V 
2{ns) x 



- COS (27lps/h 



'»}] 



p*Vh*+k*+i*r 1 (h+k+i (h-k)(h-i) 



w 



2Uh L«A z l h h°- 

1 hk - 2kl + Ik . , 1 kl 



COS 2«/r 



$* 



sin 2*h - —. (1 - cos 2»/0 

^ 2^ h* ^ \ 



(26) 



where i/j = (jtpsVh 2 + k 2 + l 2 )/h. It will be remembered 
that k, k, I are here the positive values of the indices 
arranged in the order h > k > I. 

The asymptotic value of I(S), that is the value about 
which it oscillates for large values of s, is given by the 
first term only of equation (25), 

V'(o) 



I(S)~- 



2(jw) 8 t/' 



(27) 



so 7(5) falls off as s~ 2 . The coefficient of s~ 2 is positive, 
since V(t) has its greatest value for t = o, and V (t) must 
be negative for t = o. In fact, it is easily seen that - V'(o) 
is the area of the projection of the crystal on the reflecting 
plane. Equation (27) may therefore be written 

I(S)~A/zU(tu)\ ■ . . (28) 

where A is this area. An alternative expression for the 
^ymptotic behaviour of I(S) may be got from equation (20). 
For s large sin (nTs) undergoes rapid alterations in value 
w ith changes in T (that is, with change in position in the 
Cr ystal), and in the integration over the volume of the 
Cr ystal sin 2 (tiTs) may be replaced by its average value, 
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Equation (20) becomes 



zU(m) 

which is an interesting contrast to equation (7). The 
value of I(S) for s small depends on the volume average 
of the thickness of the crystal measured perpendicular to 
the reflecting planes ; its value for s large depends on the 
volume average of the reciprocal of the same quantity. 

It is not to be expected that the small oscillations about 
the mean value caused by the trigonometric terms in the 
expressions for the line-profile would ever be observed. 
Even if all the crystals forming the powder specimen had 
the same shape it is unlikely that they would all have the 
same size, and variations in the size would cause the oscil- 
lations to come at different places and so smooth out the 
line-profile. 
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CHAPTER V 

MISTAKES IN LAYER STRUCTURES : GENERAL 
CALCULATION 

MISTAKES IN THREE DIMENSIONS: POWDER 
PHOTOGRAPHS 

1. So far it has been possible to treat each unit cell of 
a crystal as a single scattering centre, and all unit cells as 
identical. If, however, the unit cells are not identical, it 
is necessary to consider the actual distribution of atoms in 
the cell and the combination of scattered rays of different 
amplitudes. The effect of the distribution of atoms can 
be summed up in a quantity called the ' structure ampli- 
tude '. The amplitude of the X-rays scattered by each 
atom can be calculated as a function of the scattering angle 
if the distribution of electrons in the atom is known. It 
is more or less proportional to the atomic number, and is 
greatest at low angles (James, pp. 65-72). It is generally 
denoted by/, with a subscript to distinguish different atoms. 
Tables giving /as a function of sin 0/?. have been prepared 
for most atoms {Internationale Tabellen, vol. II, pp. 571-4 ; 
^ ra gg> J 933» PP- 3 2 9~34)- ^ f 1S known for each atom, 
the effect of their different positions within the unit cell 
can be determined. Suppose that there are atoms of scat- 
tering factors f v / 2 , "•, /„ at positions in the unit cell x lt 
yv*i] x 2 ,y 2 ,z 2 ; -; x„,y n , z n \ where x, y, z are given 
as fractions of the edges of the unit cell. For the hkl re- 
flection the rays scattered at (1, 1, 1) would be advanced 
in phase by 27t(h + k + I) radians with respect to those 
scattered at (o, o, o), so the rays scattered at (x, y, z) are 
advanced in phase by zn{hx + ky + Iz) radians. The re- 
sultant of the rays scattered by the various atoms can be 
constructed graphically as shown in figure 1. Each atom is 
•"^presented by a line of length/at an angle of 27t(hx+ky +lz) 
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to the horizontal, and the resultant F is represented in 
magnitude by the length of the line joining the beginning 
of f x to the end of /„, and in phase d by the angle that 
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!hxa*k»t*k*) 



O A 

Fig. 1. — Graphical construction for the structure amplitude. 

this line makes with the horizontal. Trigonometrically 
expressed, 

F°- = A 2 + B\ . . . (1) 
cos 6 = A/ \F\, sin 6 = Bf \F\ t . . (2) 
where 

A =f l cos27i(hx l +ky 1 +lz 1 ) + / 2 cos 27t(hx 2 +ky 2 +lz 2 ) 

+ ••■+/„ cos 2n(hx H + ky n + lz„), . . (3) 

B =f 1 s'm27i(fix 1 +ky 1 +lz 1 ) + / 2 sin 2Ji(Iix 2 +ky 2 +lz 2 ) 

+ •••+/„ sin 27i{hx n + ky H + lz n ) . (4) 

These expressions may be simplified if the positions of some 
or all of the atoms arc related by the symmetry elements of 
the space group. For example, in many cases the origin 
of the unit cell is a centre of symmetry, and an atom at 
the point (x, y, z) implies another at the point (1 - x, 1 - w 
1 - z). In the equation for B, therefore, for every term 
/sin 2n(hx + ky + Iz) there is another 

f$m27t{fi{i-x)+k(i-y)+l(j-z)} = -/sin 27i(hx+ky+ls), 



which cancels it, so that B is zero and d is a multiple of n. 
In such cases it is not really necessary to consider the phase 
angle. Since the lines representing the structure ampli- 
tudes are all horizontal, differing merely in direction, they 
may be added algebraically instead of trigonometrically. 
They are taken as positive if S = o, negative if <5 = n. In 
these cases then F = A, B = o. Even if B does not vanish 
for all reflections it is always possible to choose a new origin 
in the unit cell so that it vanishes for the reflection in which 
one is immediately interested. A shift of origin to the 
point (£ , f] } Q would decrease all the angles in figure 1 by 
27i(h£ + kr) + /C) ; that is, it would merely rotate the figure 
through this angle. Such a rotation would not change the 
magnitude of /'', but the angle d could be reduced to zero 
by a proper choice of the new origin. 

The simplifications produced by all the possible com- 
binations of symmetry elements have been tabulated (Inter- 
nationale Tabellen, vol. I ; Lonsdale, 1936). It may 
happen that for some values of h, k, I both A and B are 
zero, and hence that there is no diffracted beam even though 
the conditions for reflection derived in chapter II are satis- 
fied. The corresponding point of the reciprocal lattice is 
then unaccompanied by any region of X-ray diffraction. 
Such reflections of zero intensity are called absences. They 
are of great help in fixing the symmetry of the crystal, but 
care has to be taken not to confuse absences due to the 
symmetry of the crystal with ' accidental ' absences, which 
occur when F 2 for a particular reflection happens to be too 
small to be observed. We are, however, concerned not 
with the determination of structures, but with the reciprocal 
space corresponding to a crystal in which F has different 
values in different unit cells. 

2. Mistakes are particularly frequent among the ' layer 
structures '. In these the crystal can be regarded as con- 
sisting of a number of layers, the forces holding each layer 
together being usually much greater than those holding the 
separate layers together. Examples are some silicates such 
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as the micas, cadmium iodide, graphite, and hexagonal 
cobalt. If the forces between the layers are weak the layers 
may be displaced parallel to themselves, the arrangement 
of atoms within the layer being maintained, but the proper 
sequence of atoms in passing from one layer to the next 
being impaired. Even if the forces between the layers are 
strong, there may be several positions in which one layer 
can fit on to that below it without much difference in energy, 
and so occasional mistakes are possible. Cobalt is an 
example of this. It consists of close-packed planes of 
atoms, such as those shown in figure 2. If a plane of atoms 

o o o o I 

0000 

°o®d£o°o • 1 

o u ® u e^e o I 

Fig. 2. — Possible positions of planes of atoms in close packing. 

has the position A, the succeeding one can take either the 
position B or the position C. Suppose it takes the posi- 
tion B. The third plane could take either the position C 
or A. If the cobalt is in the hexagonal modification it 
should take the position A, and successive planes should 
maintain the sequence ABABABABA---. As far as the 
atoms in contact are concerned, however, the combination 
AB can be followed equally well by A or C, and it would 
not be surprising if an occasional plane forgot to consult 
the atoms with which it was not in contact, and settled 
down in the wrong place. Photographs of hexagonal cobalt 
do show peculiarities most easily explicable by mistakes of 
this kind. The layer silicates have much more complicated 
structures, but many of them also have three (or perhaps 
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six) possible positions of equilibrium. With some of them 
there is another possibility of mistakes. Layers of water 
or other substances can be adsorbed in the spaces between 
the layers, thus increasing the interlayer spacing. If the 
number of layers of water molecules is the same in every 
interlayer space the crystal remains perfect, but if the num- 
ber of layers of water molecules varies there will be erratic 
variations in the interlayer spacing. 

3. Let us first consider a hypothetical case in which the 
a and b axes lie in the plane of the layer and c is perpen- 
dicular to it, and the mistake consists of a shift of the layer 
by half the b axis. If F is the structure factor of a cell in 
an unshifted layer, the structure factor of a cell in a shifted 
layer is easily seen to be ( -) k F, since the ^-coordinate of 
every atom is increased by 1/2, and both A and B are 
changed in sign for k odd, and are unchanged for k even. 
If we try to calculate the amplitude of reflection by the 
method of chapter III, paragraph 4, we find that the sum- 
mations for lines parallel to the a and b axes in the n th 
layer go through successfully, leading to the expression 



„ sin (nXu/a) sin (nYv/b) 



71U 



nv 



where F n is F or ( -) k F depending on whether or not the 
«"' layer is a shifted layer, but the third step, of adding 
the amplitudes diffracted by successive layers to give the 
dependence on w, cannot be taken, since the layers are not 
identical. As we are considering the effect of mistakes on 
reciprocal space, not the effect of small particle size, we 
may suppose that X t Y and Z are large, and the result 
therefore means that the reflection is practically concen- 
trated into a region with very little extension in the direc- 
tions of a* and b*, but that there may be a considerable 
spread in the direction of c*. We have now to find some 
Way of estimating this. 

Suppose that the crystal starts with a layer in the ' right ' 
Position, and that each added layer tends to be in the right 
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position, but that there is a chance a of its going into the 
' wrpng ' position. If it goes into the wrong position, 
every layer after it will go into the wrong position, until 
there is another mistake, which restores the layers to the 
right position. Suppose that R„ is the probability of 
the n th added layer being in the right place, and W n the 
probability of its being in the wrong place. Clearly 

R„ + W H =i . . . ( S ) 

The probability of the (n + i) th added layer being right 
can be calculated as follows. If the n xh added layer is 
right the (n + i) ,h will be right unless there is a mistake, 
and if the n lh added layer is wrong the (n + i) ,h will be 
right only if there is a mistake. The probability of the 
(n + i) th being right is therefore the probability of the n* 
being right multiplied by the probability of no mistake, 
plus the probability of the « ,h being wrong multiplied by 
the probability of there being a mistake. In symbols this 
becomes 

R„ +1 = (i - a)R n + *W H 

= (i - a )P„+a(i -P„) 

= a+(i -2a)P n ... (6) 

This equation may be solved by trying a solution of the 
form 

Rt-P + Qp* ... (7) 

where P, Q and p arc constants to be determined. Sub- 
stitution gives P = !,/> = i -2a, so 

R n = * + 0(1 - 2a)» . . (8) 

Since R is i, the starting layer being correct, Q must 
have the value J, and 

R n = |{i + (i - 2a)"} . . (9) 

In the next paragraph the average value of the product of 
the P's of two layers separated by n interlayer spaces will 
be needed. This is clearly P 2 multiplied by the proba- 
bility of the two layers being the same (that is, both right 
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or both wrong), plus (-)*P 2 multiplied by the probability 
of the two layers being different. A little consideration 
shows that these probabilities are R n and W n , so that the 
required average is 

P 2 P H+ (-)*P 2 (i -R n ) 
= P 2 , k even, .... (10) 
= (1 - 2a)"P 2 , k odd. . . (11) 

4. Let us now try to find the intensity of reflection for 
w = o, that is at the point where the intensity is concen- 
trated for a perfect crystal. The amplitude is given by 

G =P 1+ P 2 +P 3 + ••• + P V , . . (12) 

where N = Z/c and P,„ the structure factor of the n th 
layer, is equal to F for k even, and is either F or - F for 
k odd. If k is even therefore G = NF = ZF/c, just as for 
a perfect crystal. If k is odd the terms in equation (12) 
will be positive for a while, then negative for a while, 
positive again, and so on, so that the sum is indeterminate. 
It is, however, possible to find its probable value, or rather 
its probable value without regard to sign. Depending on 
how the mistakes occur, G could have any value between 
- NF and + NF, and its average value will be zero, but 
usually it will have some finite sum either positive or nega- 
tive. In calculating the X-ray intensity the sign does not 
matter, and it is the average value of G- that is wanted. 
Squaring equation (12) gives 

G 2 = F* + Fi + F* + - + P v 2 

+ 2(P 1 P 2 + F 2 F 3 + F 3 F, + ». + Fjv-xP.v) 

+ 2(*iF, + F 2 F, + F S F B + .» + Pj^A) . 

+ ••• 

+ 2F l F N (13) 

where the terms have been arranged so that in the first 
'uie the two P's multiplied together refer to the same layer, 
'n the second they belong to adjacent layers, in the third 
there is one intervening layer, and so on. Suppose that 
Jm is the average value of the product of the P's of two 
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layers separated by m interlayer spaces. . Then 

G 2 = NJ Q + 2{N - l)/, + 2(N - z)J 2 + - + */*_! . (14) 

For the particular mistake considered in the last paragraph' 
J m was found to be 

/„-(! -2am • . • ("J 

so 

G 2 = F 2 [ 2 {N + (N - i)(i - 2a) + (TV - 2 )(i - 2a) 2 

+ - + (1 - 2a)*- 1 }-N]. . (15) 

The expression within the curly brackets is a combined 
arithmetic-geometric series, and it is not difficult to find 
its sum. In the present case, however, N is large and 
(1 - 2a) m decreases rapidly, so it may for practical pur- 
poses be regarded as geometric series of the form N(i — 2a)* 
summed to infinity, and 

G 2 = 2 NF 2 {i + (1 - 2a) + (1 - 2a) 2 + -} 

= 2NF 2 /2«. 

= F 2 Z/co> (16) 

The total intensity of the reflection is F 2 Z/c, so its integral 
breadth in w is a. Thus the greater the chance of a mis- 
take the more the region of appreciable intensity of reflec- 
tion is spread out along the c* axis in reciprocal space. 

5. Having found the maximum intensity of reflection as 
a function of position in reciprocal space, that is G 2 (w) =I(zo) 
for w = o, let us try to evaluate I(w) for an arbitrary value 
of w. In this case each layer will be reflecting 2nw out of 
phase with the preceding one, and the intensity of reflection 
as a function of w will be given by 

I(w) = G 2 (w) = A 2 + B\ 
where 

A = F t + F 2 cos 27iw + F 3 cos qjvio + F 4 cos 671a; + •••, I 
B = F 2 sin 27ZW + F 3 sin <\nw + F 4 sin (mw + •••. 

By squaring and arranging terms as in equation (13) above 
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this becomes 
I(tv)=F 2 +Fl + FI + - 

+ 2[F X F 2 COS 27ZW + F 2 F 3 (C0S 27CZ0 cos 4?rw 

+ sin 27iw sin ^mo) + •••] 
+ 2[-F 1 F 3 cos 4JK0 + F 2 F 4 (cos 2tvw cos 6tw 

+ sin 27iw sin 6jra>) + •••] 
+ ••• 
= F 2 + F\ + Ff + - 

+ 2 cos 27tw(F 1 F 2 + F 2 F 3 + F 3 F 4 + •••) 
+ 2 cos ¥vw{F x F z + F 2 F A + F 3 F 5 + •••) 
+ 2 cos 6jto;(F 1 F 4 + F 2 F S + F 3 F 9 + — ) 

+ - • • (17) 

On replacing the products of the F's by their mean values 
this becomes 

I(w) = NJ + 2(N - l)/j COS 271W + 2{N - 2)J 2 COS \JVU1 

+ 2{N - 3)/ 3 cos 6mo + — , . (18) 

where the f's have the same meaning as before. This series 
also could be summed strictly for J m = F 2 (i - 2a) m , but 
for the present purpose it suffices to approximate the sum 
by an integral : 

/•CO 

I(w) '— 2N1 J in cos 27iwm dm . . . (19) 

Jo 

= 2iVF 2 l (1 - 2a)'" cos 27twm dm . (20) 
Jo 

= 2iVF 2 | e- ym cos 27iwm dm . . (21) 
Jo 

where y = - log (1 - 2a) =. 2a. 

This integration may be done by parts : 

1 



e~ ax sin bx dx 



r 1 1°° bc K 

I e-" x cos bxdx = — e - "* cos bx\ — I 1 
Jo a J aj 

ib "I °° 6 2 f c0 

= + — e~" r s\nbx\ - I e' ax cos bxdx, 
a a- J a 2 ) 



54 
so 



X-RAY OPTICS 



Ie-™ cos bx dx = -, 
o « 



f 



e~ ax cos bx dx = a/(a 2 + b 2 ). . (22) 



In the present case a =y,b = zttzv, so 

/(.) = - 2/ ^v . . 

W y 2 + (27TW) 2 

. = aZF 2 
" c(a 2 + 7i 2 «! 2 ) 
For to = this reduces to equation (16), as it should, 
course of I(zo) as a function of zo is shown in figure 
is very similar to that of sin 2 7iNu/(tiu) 2 , 
shown in figure 2 of chapter III, except • ^ • 
that it does not oscillate. Both go to 



(23) 
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Fig. 3. — Variation of intensity as a 
function of zv for a layer structure 
with mistakes. 



Fig. 4. — Portion 
of the recipro- 
cal lattice for 
a layer struc- 
ture with mis- 
takes. 



zero as the inverse square of the variable for large values 
of the variable. The broadening of the reflections differ* 
from that due to a small thickness of the crystal in the 
c direction in that only those reflections with k odd arc 
affected. A section of the reciprocal lattice is shown il 
figure 4. 



It is perhaps worth while to check that the total intensity 
of the reflection is actually ZF-/c, as this gives some assur- 
ance that the approximations used have not introduced 
serious errors. The total intensity is 

/ = f I(w) dw = (v.ZF 2 /c) f + (a 2 + nW) l dio 

= (a.ZF 2 /c)(no)- 1 arc tan (tiw/u.) . 

At the upper limit nw/a has a large positive value, so that 
arc tan (tcw/o.) is practically 71/2, and at the lower limit it 
has a large negative value, so that arc tan {nw/a.) is prac- 
tically - 71/2. The total value of the integral is therefore 
ZF 2 /c. 

6. In the foregoing paragraphs the calculation of the 
intensity of reflec- 
tion from a layer 
structure containing 
mistakes has been 
carried through in 
detail for a particu- 
larly simple variety 
of mistake. The 
method adopted is, 
however, of general 
application in the 
calculation of the 
^tensity of reflec- 
tion from a layer 
structure containing mistakes. If A v A 2 , •••, B v B 2 , •••, 
ar e the A and B parts of the structure factors of the unit 
Ce lls in the various planes (equations (i)-(4) above), and c* 

ls taken perpendicular to the layers, the contribution of 
l he M th i ayer t0 the £ part £ the tQtal am pii tUf j e f 

re flection is (figure 5) 

A n cos 27t(n - i)w - B„ sin 2n(n - i)w, 




Fig. 5. — Construction for the A and B 
parts of the amplitude reflected by 
the w lh layer. In this figure ■zmnu 
should be 27t(h — i)iv. 
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and the B part is 

A n sin 27i(n - i)w + B n cos 27i{n — i)zo, 

and the intensity of reflection as a function of w is, the 
terms being arranged as before, 

I{w) = G\w) = A-{w) + B°-(w) 
= (A x + A 2 cos27iw + ^ 3 cos4.tw + y1 4 cos6?iw + ••• 

- B.> sin 27T.W - B 3 sin 471W - B i sin 6tcw - •••) 2 
+ (A 2 sin 27iw + A 3 sin j\ttzo + A A sin (mzv + ••• 

+ B x + B 2 COS27IZ0 + B z COSi\7lW + B a cqs6tiw + •••)§ 

+ 2(A V A 2 + B,fi 2 ) COS 2-TO 

+ 2(^2-^3 + fi 2 5 3 )(cos 27rw cos 47EW + sin 271ZV sin 4rra>) 

+ 2(A 2 A i + B 3 B 4 )(cos 4.TW cos 6ttw + sin 4^20 sin 6rrw) 

+ ••• 

+ 2(A 1 A 3 + B X B 3 ) cos 4jrzt' 

+ 2(A 2 A 4 + B 2 B .,)(cos 2^w cos 6rra.> + sin 27110 sin 6rra') 

+ 2(A 3 A b + 2? 3 B 5 )(cos 4^rw cos 8tiw + sin jjiw sin 8;rw) 

+ ••• 

+ ••• 

- 2(A l B 2 - A 2 B X ) sin 27izo 

- 2(A 2 B 3 - ^ 3 fi 2 )(sin qrno cos 2nw - cos %mo sin 2706) 

- 2(A 3 B A - A 4 B a )(s\n dnw cos tpno - cos 6tcw sin 4700) 

- 2(A l B 3 - A 3 B X ) sin 4;™ 

- 2(A 2 B A - ^ 4 R 2 )(sin 6tt.w cos 27iw - cos 6tko sin 27110) 

- 2(A 3 B 5 - ^ 3 B 3 )(sin Stoo cos 47ml - cos Stzw sin 4.-774') 



= {A\ + Bf) + {A 1+ ££) + {A'i + B%) + - 
+ 2 cos 27iw[{AiA 2 + B X B 2 ) + (A 2 A 3 + B 2 B 3 ) 

+ {A 3 A t + B 3 B t ) + •••] 
+ 2 cos 47ro)[(^ 1 ^ 3 + B y B 3 ) + (A 2 A 4 + B 2 B A ) 

+ (A 3 A 5 + B 3 B 5 ) + •••] 
+ 2 cos brno^A^ + 5^4) + (^2^ s + 5265) 

+ (^ 3 ^ 6 +5 3 S 6 ) + •••] 
+ ••• 
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- 2 sin 27iw{{A x B 2 - A 2 B y ) + (A 2 B 3 - A 3 B 2 ) 
+ (A 3 B i -A i B 3 ) + -] 

- 2 sin 4,-rop^ - A 3 B X ) + (A 2 B A - A 4 B 2 ) 

+ (A 3 B 5 -A 5 B 3 ) + -] 

-- (25) 

This equation contains terms of two kinds, those of the 
form A A' + BB', and those of the form AB' - A'B, where 
A and B belong to one layer and A' and B' to another 
layer. If J m is the mean value of terms of the first kind 
when the layers are separated by m interlayer spaces, and 
K lh the mean value of terms of the second kind, equation 
(25) becomes 

I{w) = NJ +2(N-l)J l COS 271Z0+2(N-2)J 2 COS ^7CW\ 

+ - + 2/.v_, COS 27C(N - l)tV ^ , fi v 

-2(N-i)K 1 sin 27iw -2(N -2)K 2 sin <\mo | 
_ ... _ 2K y _i sin 27i(N-i)w. 

This is the generalized form of equation (18), and is valid 
in every case when the number of layers is great enough 
to allow averaging the values of A A' + BB' and AB' - A'B. 
If in addition the chance of a mistake is large compared 
with i/N, but small compared with 1/2, the sums may be 
approximated by integrals, so that 

(N - m)J m cos 27iwm dm 

r- v 

- 2 1 (N - m)K m sin 27110m dm. . (27) 
Jo 

The correct upper limit is N = Z/c, but with the chance 
of a mistake within the range stated it can just as well be 
niade infinity, and (N — m) can be replaced by N. In 
niany cases therefore equation (26) may be replaced by 



/(«;) = 2Nfv n 
Jo 



cos 27iwm - K m sin 27110m) dm, (28) 



but it must be remembered that this is an approximation, 
"^d if the chance of a mistake is very large or very small 
e quation (26) must be used as it stands, or approximated 
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in a way appropriate to the particular problem. From the 
definitions of J m and K m it is obvious that J_,„ = J mt and 
K_ m = - iC, n , so that equation (28) may be written more 
symmetrically as 



m 



r +7 



(J m cos znwm - K m sin znwm) dm 



or 



/(«,) = 



(30) 



/(w) = Nj(J m cos 2ot»i - /C B , sin 27vwm) dm . (29) 
for short. In very many cases, as in the example used in 
paragraphs 4 and 5, K m is zero and the terms involving it 
vanish. In most of the problems to be considered in the 
following chapter J m has the form J e~ y "\ as in equation 
(21) above, and equation (28) becomes 

y- + {ztiw) 2 

The maximum value of this is 2NJ /y, and the total 
intensity is NJ Q , so that the integral breadth in w is y/a.1 
7. The general effect of mistakes in layer structures will 
now be obvious. Reflections whose structure amplitude 
is not affected by the mistake in question (for example, all 
reflections with k even in the example considered above) 
will be unaffected, and the spots or lines on the X-ray 
photograph will be as sharp as if the crystal were perfect. 
On the other hand, reflections whose structure amplitude 
is affected by the mistake will be drawn out in reciprocal 
space in a direction perpendicular to the layers. The 
amount of the spread will depend on the frequency of the 
mistakes and the amount of the change in the structure 
amplitude ; in the simple example considered in para- 
graphs 4 and 5 all reflections with k odd were drawn out 
into streaks in reciprocal space of length of the order of 
c*/ol. The corresponding spots on the X-ray photograph 
will in general be drawn out also, or the lines in the powder 
photograph broadened, but the magnitude of the effect 
depends on the angle that the streak makes with the 
surface of the sphere of reflection, and is most readily 
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discussed in connection with particular examples, as is done 
in the following chapter. If K m is not zero the reflection 
is displaced as well as broadened. This is readily seen 
from any of equations (26) to (29), since /(- w) is not equal 
to I{w) unless the terms involving sines vanish. One or 
two examples of displacement arc discussed in chapter VII. 
8. Besides layer structures with mistakes, which may be 
described as mistakes in one dimension, examples are known 
of substances exhibiting mistakes in three dimensions. One 
of these is an alloy of gold and copper, in the proportion of 
one gold atom to three 
copper atoms. The nature 
of the mistakes may be 
understood from figure 6, 
which represents one unit 
cell of the crystal. The 
cell is cubic and contains 
four atoms, distinguished 
by the letters A, B, C, D. 
(At first sight there appear 
to be twelve atoms, but 
only one-eighth of each 

atom marked D is really inside the cell, the remaining 
seven-eighths being in neighbouring cells. Similarly, 
only one-half of each atom marked A, B, or C is inside 
the cell.) On the average, one of these four atoms will 
be gold and three will be copper. At high temperatures 
the gold and copper atoms are mixed at random, and 
powder photographs show only the reflections corre- 
sponding to a face-centred cubic structure. At lower 
temperatures the gold atoms, behaving as if they were 
larger than the copper atoms, avoid one another, and in 
a n ideal crystal would all settle into position A, or all into 
Position B, or all into position C, or all into position D, 
the three positions not occupied by gold atoms being 
occupied by copper atoms. Thus at high temperatures 
^y position, say A, is occupied on the average by one- 




FlG. 6. — Unit cell of AuCu 3 . 



6o 



X-RAY OPTICS 



quarter of a gold atom and three-quarters of a copper 
atom, and at low temperatures either always by gold or 
always by copper. There are therefore at low temperatures 
four distinct types of unit cell : 

A, with gold atoms in the centre of the face perpendicular 

to the a axis, 

B, with gold atoms in the centre of the face perpendicular 

to the b axis, 

C, with gold atoms in the centre of the face perpendicular 

to the c axis, and 

D, with gold atoms at the cube corners. 

When the gold atoms have sorted themselves out in one 
of these ways throughout a crystal it is said to be ' ordered ', 
and a large number of extra X-ray reflections appear, cor- 
responding to the decrease in symmetry from statistically 
face-centred cubic to simple cubic. If the ordering is 
complete, each gold atom is entirely separated from con- 
tact with others by a screen of copper atoms. A com- 
paratively small increase in temperature is able to overcome 
the gold atoms' dislike of one another's company, so that it 
would not be surprising if within an ordered single crystal 
mistakes occurred such that the type of order varied from 
place to place. This is found to be the case, and by proper 
heat treatment it is easy to put the alloy into a state in 
which the reflections common to the ordered and disordered 
state are sharp, but those characteristics of the ordered state 
only are more or less diffuse. One or two other substances 
exhibiting mistakes of this type are known, but they are 
much rarer than layer structures with mistakes, perhaps 
because it is only possible to build up three-dimensional 
mistakes with exactly four types of unit cell (Bragg, 1940). 
9. The intensity of reflection from such a structure as 
a function of position in reciprocal space could be treated 
by methods similar to those of paragraph 6 above, but the 
double series there, already sufficiently clumsy, are replaced 
by sixfold series, and become almost unmanageable. Con- 
sideration of this problem will therefore be postponed to 
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chapter VII, where more advanced mathematical methods 
are introduced. Line broadening in powder photographs 
may, however, be treated here by essentially the same 
argument as was applied previously in connection with 
particle size broadening (chapter IV). To calculate the 
integral breadth of the hkl reflection, we calculate the in- 
tensity scattered with S between S and S + dS by a single 
line of unit cells parallel to the radius vector joining the 
point h, k, I of the reciprocal lattice to the origin, and sum 
for all parallel lines. If paragraph 6 is re-read, it will be 
seen that the calculation of the intensity scattered by a 
single line of cells goes through in exactly the same way 
as the calculation of I(zo), except that J m and K„, axe the 
mean values of A A' + BB' and AB' - A'B for cells separ- 
ated by m repeat distances parallel to the line, instead of 
perpendicular to the layers. With axes chosen as in 
paragraph 3, chapter IV, one may rewrite equation (27) as 

rT 
//(«) du = (2/a) \ (T - t)J(t) cos {zntu/a) dt du 

Jo 

- (2/a) [ {T - t)K(t) sin {zntu/a) dt du . (31) 
Jo 
where a is the repeat distance in the direction of the line, 
T = Na is the length of the line, t = ma, and u is written 
in place of w to conform to the notation of the previous 
chapter. The interpretation of J(t) and K(t) is fairly 
obvious ; they are the mean values of AA' + BB' and 
AB' - A'B for cells separated by a distance t along the 
line. Equation (31) is analogous to equation (3) of 
chapter IV, and I(S) is obtained from it in the same way : 

I(S) = {a/U)^H(u)dydz 

= (2/I/) [ f f (T - t)J(t) cos (27itu/a) dt dy dz 
- (2 /U) [ f [ T (T - t)K(t) sin (zntu/a) dt dy dz . (32) 
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By equation (15), chapter IV, the volume common to a 
crystal and its ' ghost ' shifted a distance t in the hkl 
direction is 

V(t)=SKT-t)dydz, 



so 



7(5) = (a/CT) [ + V(t)J(t) cos {zntu/a) dt 

- (z/U) { + V(t)K(t) sin {zntu/a) dt 

= (z/'U) fV(/X/(0 COS Znst dt 

-(z/U)[ + V(t)K(t)sinzmtdt . (33) 
Jo 
where the integration extends to a large enough value of t 
for V(t) to vanish, and s = u/a is the actual distance by 
which 5 misses the point //, k, I of the reciprocal lattice. 
It will be remembered that u is the fraction of the distance 
i/a between the reciprocal lattice points by which S misses 
the point of the reciprocal lattice. Since J{t), V(t) and 
cos znst are unaffected by a change in the sign of t, and 
K(t) and sin znst both change sign, equation (33) may be 
written 

I{S) = U-^V{t)J{t) cos znst dt 

- } V(t)K(t) sin znst dt] . (34) 

where the limits of integration are the negative and positive 
values of t for which V(t) vanishes. This equation is very 
similar to equation (21) of chapter IV, and in fact reduces 
to it if there are no mistakes, for then J(t) is constant and 
K(t) is zero. 

Ordinarily J(t) will have the form 

7(*)-/o«p(-yM)> 

where y is a function of h, k and /, and K(t) will be zero- 
If mistakes are sufficiently frequent to cause line broaden- 
ing J(t) will become inappreciable before V(t) departs 
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greatly from its value for t = o. Equation (34) then 
becomes 

exp(-y|/| ) cos znst dt 

—CO 



JoV 



zy 



■ (35) 



U y 2 + (zns) 2 
The maximum value of this is I H = zJ V/Uy, so that the 
apparent particle size is 

e =/ s // = 2/y (36) 

and the integral breadth is 

= If e cos m Xy/z cos 0. . . (37) 
When K(t) is not zero one may expect displacement as 
well as broadening of the lines. 



REFERENCES 

W. L. Bragg, 1933. The Crystalline State. London : Bell and 

Sons. 

1940. Proc. Phys. Sec., 52, 105. 

S. B. Hendricks and E. Teller, 1942. J. Client . Phys., 10, 147. 
Internationale Tabellen zur Bestimmung von Kristallstntkturen, 

1935. Berlin : Borntrager. 
L. Landau, 1937. Phys. Z. d. Sotc. Un., 12, 579. 
1. M. Lifschitz, 1937. Phys. Zeit. d. Sow. Un., 12, 623. 

1939. J. Exp. Theor. Phys. USSR., 500. 

K. Lonsdale, 1936. Simplified Structure Factor and Electron 

Density Formulae for the 230 Space Groups of Mathematical 

Crystallography. London : Bell and Sons. 
A. J. C. Wilson, 1942. Proc. Roy. Soc, A, 180, 277. 
1943- Proc. Roy. Soc, A, 181, 360. 

(References to experimental work on substances mentioned 
above will be found in the next chapter.) 



CHAPTER VI 

MISTAKES IN LAYER STRUCTURES : SILICATE! 
COBALT, GRAPHITE 



MISTAKES IN THREE DIMENSIONS: 
PHOTOGRAPHS OF AuCu 3 



POWDER 



- 



i. The mistake specially considered in paragraph 3 
the preceding chapter does not seem to occur, but a 
similar one with a shift of one-third or two-thirds of 
the b axis is found in some layer silicates (Hendricks, 
1940). Suppose that the three types of layer are A, B t 
and C, and that the crystal starts off AAA-", but that 
at each added layer there is a chance a of a mistake 
occurring, and the sequence continuing BBB"' or CCC"*i 
Suppose, further, that either mistake is equally likely. 
Let the probability of the n th added layer being A be 
A nt of being B be B n , and of being C be C B . 

Then A n+1 = (1 - a)A n + h<xB n + JaC n , . (1) 

the factors of one-half coming in because, if the n lh added 
layer is B or C and a mistake occurs, the (n + i) th will 
be A only half the time. Since A n + B n + C n = 1, 
equation 1 may be written 

A n +i = (1 - S*)4i + l*- • • ( 2 ) 
This may be solved in the same way as equation 6 of 
the preceding chapter, the solution being 

^.-l + flfc-W • • (3) 

where Q is an arbitrary constant. Similar equations can 
be set up for B n and C n , with the same solution. Choosing 
the constants so that A = 1, B = C = o, gives 

A n - *[i + 2(1 - fa)"],] 
B„ = \[x-{i-Un \ . . 
C n = |[i - (1 - fa)-]. J 
If the structure amplitude of the A plane is F, that for tto 
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, 



(5) 
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5 plane will have an A part of F cos (271k/ 2) and a B part 
of F sin (271k/ 2), and that for the C plane will have an 
A part of F cos (47^/3) and a B part of F sin (47^/3). 
Finding the mean value of the product AA' + BB' is 
simplified by the fact that B = 0, so that J m is the mean 
value of A A', and is given by 

/,„ = F*[A m + B m cos (xrck/3) + C m cos (47^/3)] 
= iF 2 [{x + cos (27lk/2) + cos (4^/3)} 

+ (1 - fa)"' {2 - cos (27ik/2) - cos (47^/3)}] 
= F 2 , k = 3 x an integer, 

= F 2 (i - 2<x)'", k = (3 x an integer) ± 1. 

In deriving equation (5) it was tacitly assumed that one of 
the planes concerned was an A plane, but the same result 
is obtained if a B plane or a C plane is assumed. It is 
readily seen that K m is zero : 

Kn = F 2 [A m x o + B m sin (2.^/3) + C m sin (4^/3)] 
= 0, 

since B m = C w and sin (47^/3) = - sin (271^/3). 

Since J m = F 2 the reflections with k/2 integral are 
unaffected, but those with k/2 fractional are drawn out 
into streaks in reciprocal space parallel to the c* axis, the 
variation of intensity along the streaks being given by 
equation (23), chapter V : 

2NF 2 y 



/(«) 



y 2 + (27iw) 2 ' 
I2NF 2 0L 

9a 2 + (pizv) 2 



(6) 



since y = 30C/2. On a rotating-crystal photograph these 
will show as streaks along curves in the photograph corre- 
sponding to lines parallel to the c* axis, but on a powder 
Photograph they will be somewhat fore-shortened, and 
hnes for which OQ is perpendicular to c* will be practically 
u nbroadencd. For example, in an orthogonal crystal the 
''nes with / = o will be sharp, as well as those with k 
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divisible by three. A more extended discussion of the fore, 
shortening is given below, in connection with line-broad 
ing in photographs of cobalt, and of the lines with / 
in connection with graphite. Hendricks (1940) has made 
an extensive investigation of silicates showing mistakes of 
this type, but has not attempted a quantitative comparison 
with theory. 

The other mistake possible with the silicate minerals is 
an increase in the interlayer spacing at random intervals. 
If the mistake consists of an increase of the spacing by 
a fraction g, each mistake will cause an increase of the phase 
of all succeeding layers of znlg. The mean value of the 
product A A' + BB' for layers separated by m interlayer 
spaces will therefore be 
J m = F 2 multiplied by the chance of no mistake in m layers 

+ F 2 cos 27ilg times the chance of one mistake in 
m layers 

+ F 2 cos ^Tilg times the chance of two mistakes i 
m layers 

+ ••• 

+ F 2 cos zmnlg times the chance of m mistakes in 
m layers 



" 



= F 2 |(i - a)'" + ma(i - a)'"" 1 cos znlg 

+ — — — a 2 (i-a)'" -2 cos 4?ilg + ••• + a'" cos znmlg . 

In the problems so far considered K m , the mean value of 
AB' - A'B, has been zero, a considerable simplification 
being thus produced. In the present case, however, K m 
does not vanish, but is given by an expression like that fof 
/„„ except that the cosines are replaced by sines, so that 

= F 2 \ my.{i - a)'"- 1 sin znlg 

t- - ■ — j a 2 (i - a)'"- 2 sin tpilg + — + a'" sin znmlg \ 



K 
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(The term in (1 -a)'" is missing in the expression for 
K„ t because it is multiplied by the sine of zero, and so 
vanishes, whereas in J m it is multiplied by the cosine of 
zero, which is unity.) With expressions as complicated as 
these for /„, and K m it is not possible to evaluate equation 
(28), chapter V, simply by elementary means, and recourse 
must be had to methods involving complex quantities. It 
may be said, however, that K m different from zero means 
that w = o is no longer the position of maximum intensity, 
so that we may expect a change in the positions of the 
reflections as well as a broadening. 

2. In hexagonal cobalt the mistake is a shift of one-third 
in the a direction and two-thirds in the b direction, or else 
twice this (paragraph 2, chapter V). If the structure 
amplitude of an A plane is F, that for a B plane has an 

(h zk\ „ (h-k\ . _ 

A part of F cos zn[ - + - 1 = F cos zn[ - - 1, and a B 



C 

h 



FsinzTt 



3 
h - k 



while that 



3 

zk 
part of F sin zn[ - + — 

\3 3 
of a C plane has an A part of F cos zni j and a B part 

of -Fsin27t( J. The problem differs from that 

just treated in that the w th layer and the (n + i) th must be 
different ; the mistake consists of the (n + z) iU being unlike 
the 71 th . Thus if the first plane is A and the second B, the 
third ought to be A, but there is a chance a that it is C. 
Let A„ be the chance that the n th added plane is like the 
first. The probability that the (n + z) th is like the (n + i) th 
»s zero, that it is like the w th is (1 - a), and that it is like 
Neither is a, so 

A n + 2 = (1 - v)A n + o x A„ +1 + a(i - A n - A nl .i), 
or A lll2 + <iA H+1 - (1 - z*)A n = a. . (7) 

This is a nastier equation than either (6), chapter V, or 
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(2), above, but it can be solved in the same way. Let 

A n = P + Q p * t 
P + Qf** + oP + olQ p »+i - (1 - 2a )P - (1 - 2a)0p» = «, 

P = */3. 

p n+2 + ap" +1 - (1 - 2a)p n = o, 



P = 



_ - a ± Va 2 + 4 - 8a 






(8) 



It will be found that the positive sign leads to reflections 
that, with the ordinary hexagonal indices, have / even, 
and the negative sign to reflections that have / odd. Let 



- a + 



V 4 -8a 



+ a s 



- a - Va. - 8a + a 5 
Po = 2 






(9) 



Then 4. -* + &/>." +Q Po 11 . . (10) 

where O e and are arbitrary constants. They must be 
chosen so that A = 1, A y = o ; the values required are 
readily found to be 



M 



(12) 



a -(I +2 Po )h(p -Pe); 

Q = (I + 2p e )/ 3 (p e - Pl> ). 
For c/Z « a « 1 these reduce to 

f>, - l - f« = exp (- 3 a/ 2 ), & = 1/6/ 
p - - 1 + Ja = - exp (-a/2), O = i/ 2 ._ 

The values of B n and C,„ the probabilities that the »* 
added plane will be fi or C, are given by equation (10) 
with the constants chosen to make B = C = 0, 
5 X - C x = 1/2, since the first added plane is equally likely 
to be B or C. It is found that 

B„ = C n = | + fy,,* + R oPo » 



where 
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R e = (1 + 2 Po )/6( Pe - Po ) = - 1/12, 
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• (i3) 



i? = (1 + 2 Pe )/6( Pn - p,) -=:- 1/4. 

The mean value of AA' + BB' is then 

J m = F*A m + FW* cos *ft^\ + F2C m cos 2vr(-— \ 



= F 2 [i(i + 2C) + (g B + 2C/?,K'» 

+ (Q + 2CK„)p*] . 



where C has been written for cos zni - 



m 



■ (14) 

For A - £ 



divisible by three this reduces to F 2 ; for A - & not divisible 
by three it reduces to 

J m n *■»(*,>» + i P< r) 

or F*[| exp (- 3a»i/2) + ;J cos wwr exp (- aw/2)] (15) 

approximately. The value of K m is readily found to be 
zero. Reflections with h - k divisible by three are normal, 
since for them J m = F 2 ; the line-profiles for the other 
reflections are given by 

/(w) = 2iVl J„, cos 27twm dm 
Jo 

= 2NF*( {* exp (-30^/2) 

+ :] exp (- aw/2) cos ntTc) cos 27iwm dm 
3a a 

2 12 



= 2iVF 2 < 



m 



^Y +{27iwf 



♦fc 



*U 



feV + [ 27r (a, + J)] 2 

. . (16) 



-)* + [27z(«>-l)] 2 
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(The second integration is done by expressing the product 
of two cosines as the sum of two cosines.) This is a rather 
more complex expression than we have obtained previously, 
but the interpretation is straightforward. The points of 
the reciprocal lattice have broad weak peaks of integral 
breadth a = 3 a/4 and intensity NF 2 /^, while halfway 
between them there are sharper stronger peaks of integral 
breadth a = a/4 and intensity t,NF 2 /^. The reason for 
this apparent halving of the reciprocal lattice is that wc 
have been taking a unit cell as consisting of a single plane, 
whereas the normal unit cell of hexagonal cobalt consists 
of two planes. The reciprocal lattice, with the ordinary 

values of /, is illustrated in 

• f \ • \ figure 1. The points with 
* I * h - k divisible by three and 

I • I I • f / even are sharp, and with I 

I l c « * f odd are missing. The 

• q * I f • § points with h - k not divis- 
* t * i ible by three and / even arc 

• I I • I drawn out into long streaks 
5 ♦ } \ parallel to the c* axis ; those 

FlG. ..-Portion of the recip- with l odd are ^° dnuvn 
rocal lattice for cobalt. out, but only one-third as 

much. The integral breadths 
referred to the ordinary values of / are twice those given 
above, so that a - 301/2 for reflections with / even, and 
a = a/2 for reflections with / odd. 

In the case of cobalt the strict calculation of the intensity 
has been carried out, equation (26), chapter V, for the line- 
profile and equations (9), (11) and (13) for the p's and Q's 
being used. For a small compared with c/Z the reciprocal 
lattice for normal hexagonal close-packing is obtained, for 
a large compared with c/Z but small compared with one- 
half it is as just described, and for a~i/2 the fuzzy 
regions practically disappear. For « a little greater than 
one-half the p's become imaginary, but I(tv) can still be 
interpreted, fuzzy regions appearing in places correspond- 
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ing to the points of the reciprocal lattice for cubic close- 
packing. As a approaches unity these regions coalesce 
into sharp points, the final state being the superposition of 
the reciprocal lattices for the two orientations of the cubic 
crystal corresponding to the sequences of planes 
ABCABC-" and ACBACB-. This is as it should be, 
since an hexagonal crystal with a mistake at every plane 
(x = 1) would be cubic close-packed, but the beginnings 
AB and AC are equally likely. 

A careful experimental study of the broadening of the 
X-ray reflections from imperfect hexagonal cobalt has been 
made by Edwards and Lipson (1942), both single-crystal 
and powder photographs being used. The single-crystal 
photographs showed in a striking manner the extension 
of the regions of reflection along lines parallel to the c* 
axis, and the fact that reflections with / odd were much 
less affected than those with / even. They were, however, 
unsuitable for quantitative measurements, and for this pur- 
pose a powder photograph taken with crystal-reflected 
nickel radiation was used. (Crystal-reflected radiation 
is practically monochromatic, 
consisting almost entirely of 
one doublet, and so gives 
photographs with less back- 
ground than merely filtered 
radiation does.) The integral 
breadths were derived by the 
rnethod of Jones (1938). The 
'ntegral breadths are given in 
the second column of table I. 
They depend, however, on 
the actual value of /, as well 
*s on a and the parity (odd- 
ttess or evenness) of /. This is easily seen from figure 2. 
A ' spike ' of length in reciprocal space 6 corresponds 
for / = ^ to a range of sin proportional to B^C^ whereas 
*or / = / 2 it corresponds to the greater range proportional 




Fig. 2. — Dependence of the 
breadth of a reflection on /. 
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to B Z C 2 . The relation between /? and / is most easily 
from the equation 

sin 2 = X\h* +hk+ & 2 )/3« 2 + x 2 / 2 /^ 2 , . (17). 

which relates the angle of reflection to the indices 
reflection for an hexagonal crystal. Differentiation giv 



got This 
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2 sin cos 0— = — -. 

dl 2C- 



(18 



If a is the integral breadth in /, or w, which amounts to thi 
same thing, 

a d0 



dl c- sin 20 



- H dd A 2 



/ 



/<•• 






TABLE I 

Calculation of a, the probability of a mistake per added plane of 
atoms, in a specimen of Cobalt 



Line 








a 


a corr. 


ion 


0-0092 


0043 


0086 


0-088 


1012 


0-0320 


0089 


0059 


0073 


1013 


0-0256 


0051 


0-102 


o-iii 


2021 


00076 


0044 


O088 


0093 


2022 


0-0406 


0105 


OO7O 


0088 


1014 


0-0844 


0099 


0-066 


0-119 


2023 


00592 


0042 


0-084 


0093 



It follows that ft = o for / = o, and that for a given value 
of d it is smaller the smaller the value of /. The experi- 
mental values of a, derived from those of fi by means of 
equation 19, are given in column three of the table. I* 
will be seen that they fall into two groups, those with 
a ~ 0-045 (^ °dd), and those with cr~o-io (/ even)- 
The theory given above predicts a = 3<x/2 for the lines 
with / even, and a = a/2 for the lines with / odd. The 
value of a, the chance of a mistake per added plane, ' 8 
thus about 0-090 if calculated from the lines with / odd ( 
and about 0-067 if calculated from the lines with / even* 



This is not entirely satisfactory agreement. Though the 
value of a is not very small compared with one-half, 
errors due to the various approximations made (for 
example, the use of equations (12) instead of (9) and (10) ) 
do not seem sufficient to account for the discrepancy 
directly. A more probable explanation is limitation of 
the range of photometry. The theoretical profile of I(zv) 
is shown in figure 3. In calculating the integral breadth 
as a = J7 dw/I{6) it is tacitly assumed that it is possible 
to photometer a sufficient length on the film on each 
side of the line to get practically all the intensity. In 
the present case the region that it was possible to photo- 




Fig. 3. — Correction of the intensity for limitation of the range 
of photometry-. 

meter corresponded approximately to the range from 
to = - 0-15 to to = +0-15, so that the intensity corre- 
sponding to the shaded portion of the figure was lost. 
A correction can be calculated (Wilson, 1942), and it is 
found that lines with / even arc affected more than the 
lines with / odd. The corrected values of y. are given 
»n the fifth column of the table, and it will be seen that 
the discrepancy has practically disappeared. The mean 
°f the corrected values of a is 0-096 ± 0-008 from the odd 
hnes and 0-093 ± 0-014 from the even lines. It seems, 
therefore, that the assumption of mistakes in the structure, 
distributed at random along the hexagonal axis, explains 
satisfactorily the broadening of the lines and accounts 
6 
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quantitatively for the difference between the lines with 
different parity. 

The correction to the breadths may be calculated as 
follows. If the integral breadth in w is a, 

«»-#»■ • • («J 

a 2 + (xwy 

(equation 30, chapter V, with y replaced by 2a). If the 
range of photometry is from w = - x to w = + .v, the 
observed maximum intensity will be 

/; = /(o)-/(*) 

= I<fl 2 X 2 /{o 2 + 7l 2 X 2 ) . 

instead of I = /(o) = NJ Q /a. . 

The observed total intensity will be 
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V - I /(») dw - 2XI(X) 



«A7.r 
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arc tan (nx/a) 

L °" : 



7T.xa 



+ {my 



instead of / 



-J n 

J —or, 



w) dw - NJ Q . 



The observed integral breadth will therefore be 

, 2(7 r 



n 



[ (1 + a 2 /7t 2 x 2 ) arc tan (nx/a) - o/nx], (25 



, 



instead of a as it ought to be. For x =0-15, a' = 0-045 
(mean for / odd) a little trial and error gives a = 0-05 1» 
and for a' = 0-098 (mean for / even) it gives a = 0-140. 
The corresponding values of a are 0-102 and 0-093, m ^ r 
agreement with one another and with the values found by 
Wilson by correcting each line individually. 

3. Another substance exhibiting mistakes is graphite, 
and it has been the subject of numerous investigations in 
connection with lamp black, coke, arc carbons and other 
commercially important products. The ideal graphite 
structure consists of layers like those shown in figure 4t 
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Fig. 4. — Possible positions of graphite 
layers, slightly displaced to show 
hidden atoms. 



alternate layers taking the positions A and B. Well- 
crystallized graphite usually consists chiefly of this struc- 
ture, but contains 
also an appreciable 
fraction of another 
structure in which 
the layer sequence 
is ABCABC — 
(Lipson and 
Stokes, 1942). The 
relation between 
these structures is 
similar to that 
existing between 
hexagonal and 
cubic cobalt, but intermediate structures do not seem to 
have been reported, and the remainder of this discussion 
will be confined to disorder in the ideal form. The binding 
between adjacent layers is weak, so that they are readily 
displaced or rotated relatively to one another by any 
amount. This behaviour contrasts with that of cobalt and 
the layer silicates, which allow only a small number of dif- 
ferent displacements. In addition, the crystals are small, 
so that it is necessary to consider ' particle size ' broadening 
as well as ' mistake ' broadening. 

The structure of ordinary, non-ideal, graphite may be 
pictured as a sequence of planes ABAB"' for a certain 
distance, but at each added plane there is a chance of a 
random translation and rotation, destroying completely the 
regular phase relation between atoms separated by the dis- 
placement. This will not affect the 00/ reflections, so that 
for them the intensity will depend on w through the 
ordinary factor 



sin 2 jiNw 
I(w) - I' 2 — — —— 
' ' (too) 2 



(26) 
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where N is the number of unit cells in the crystal in the 
c direction (chapter III), but all other reflections will be 
drawn out parallel to the c* axis by an amount depending 
on a. The calculation of/,,, is complicated slightly by the 
fact that a unit cell consists of two layers. If a is not too 
large, however, it will be near enough to take the chance 
of a mistake per unit cell as 2a, and treat the crystal as 
consisting of N layers of unit cells, instead of 2N actual 
layers. Then the mean value of AA' + BB' is F 2 if no mis- 
take intervenes between the layers of unit cells concerned, 
and zero if a mistake does intervene, since there is no 
regular phase relation between layers separated by mis- 
takes. The value of J m is therefore F 2 multiplied by the 
chance of there being no mistake in m added layers of 
unit cells, so 



J m = F 2 (i - 2a)'" - F 2 exp (- 2am). 



(27) 



The intensity of reflection for reflections other than the 
oo/'s will therefore depend on w through the factor 



/(«) = 



a- + (tiw) 2 



■ (28) 



(equation (24), chapter V). For a not too large a specimen 
consisting of such crystals will give 00/ powder lines as 
sharp as the length of the crystal in the c direction permits, 
and hkl lines with breadths depending on a and / through 
a relation similar to equation (19) : 



P 



r-i 

c 2 sin 20 



• (29) 



the integral breadth in w of (28) being a. 

The reciprocal lattice just described, consisting of rela- 
tively sharp points on the c* axis, and all other reflections 
represented by streaks, suffices to treat powder photo- 
graphs, but it is not strictly correct. The regions separated 
by mistakes are not only translated with respect to one 
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another, but are also rotated by random amounts. This 
corresponds to a rotation of the reciprocal lattice about the 
c* axis, so that the reflections, other than the oo/'s, will be 
turned into cylinders with centres on the c* axis, the inten- 
sity being uniform around the circumference of the 
cylinders, but depending on w according to equation (28). 
A single-crystal photograph of graphite containing mistakes 
of the sort considered here would be expected to show, 
therefore, broad bands in directions perpendicular to the 
c* axis, rather than narrow streaks parallel to it. 

For large values of a the approximations leading to 
equations (28) and (29) are no longer valid, and further 
consideration is needed. It does not seem worth while to 
treat intermediate cases, but the final state of a shift at 
every added plane (a = 1) is of importance. The 00/ 
reflections will be, as before, relatively sharp points depend- 
ing on zv according to equation (26), but the other reflec- 
tions will be joined up into continuous lines parallel to 
c, and of intensity uniform in w except for the falling off 
of the atomic scattering factor with increasing sin 0/X. The 
crystal behaves essentially as if it were single for the 00/ 
reflections, giving points as sharp as its size allows, but as 
an assemblage of non-co-operating planes for all other 
reflections, giving lines in reciprocal space perpendicular 
to the planes (chapter II, paragraph 4). Crystals as faulty 
as this will not be large enough to give single-crystal photo- 
graphs, so that the rotation of the reciprocal lattice, pro- 
ducing cylinders instead of lines, need not be considered. 

A powder photograph of a specimen consisting of such 
crystals will show lines of two types. Lines of the first 
type correspond to the 00/ reflections, and are essentially 
similar to the broadened lines considered in chapter IV. 
The apparent particle size may be got from equation (26), 
but it is perhaps instructive to use the method of chapter 
IV. For this purpose the crystal may be considered 
a cylinder (not necessarily circular) of length Z, so that 
V{t), for a shift parallel to the c axis, is P(Z - \t\ ), where 
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P is the cross-section of the cylinder, so that the apparent 
particle size is 

e = V~ l jV(t) dt 
= (PZ)~ 1 i +Z p(Z-\t\)dt 

=z. . ~ 7 : . . . . ( 3 o) 

Thus the apparent particle size calculated from the integral 
breadth gives the true length in the c direction. For 
comparison with reflections of the second type, to which 
the concept of integral breadth is not applicable, it is inter- 
esting to work out the half-width of the 00/ reflections. 
For these the line profile is given by 

. sin 2 ?rZ.y 

1(S) '^F ■ ■ ■ (3,) 

where proportionality factors have been omitted, 
S = 2 sin 0/X is the distance from the origin in reciprocal 
space, and s is the amount by which S differs from its value 
S Q at the point 00/. The intensity has the value Z 2 for 
s= o, and falls to half this for the value of s satisfying 



sin 



nZs 



(Ttsy 



or 



s = ± O'442/Z 



(32) 



approximately. Thus the half breadth in s is o-88/Z, and 
the half breadth in 20, the quantity corresponding to 
integral breadth, is A/cos times that in .?, or 

0-885* 



The apparent 
breadth is 



particle size 



ZcosO 
calculated 



e' =?./fi cos = I-I3Z, 
and is thus 13 per cent large. 



• • (33) 
from the half 

• • (34) 



Lines of the second type consist of an amalgamation of 
all hkl reflections having the same values of // and k, and 
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differ from any so far encountered. The calculation of the 
line profiles involves mathematical functions beyond trie 
scope of this book, so we must content ourselves with 
a qualitative treatment. As we have seen, through every 
lattice point (h, k) in the plane / = o there passes a rod- 
like region of reflection. The intensity on the film with 
between and + dO will depend on the amount of this 
rod that lies between concentric spheres of radii S and 
S + dS, where dS = 2 cos dO/l (figure 5). The radius 
of the rod is of the order of 1/2/), where p is the square 
root of the area of a layer. Let s be the amount by which 




Fig. 5. — Calculation of the line profile for hid reflections. 

S differs from its value S at the point hko. The intensity 
due to the rod will be zero until 5 is approximately 
5 - i/2p, when it will increase rapidly until S = S Q . 
For S > S the volume of the rod intercepted by the spheres 
continues to increase until S is approximately S + i/2p. 
Beyond this point it decreases again, but less rapidly than 
it initially increased. For S > S the intensity on the film 
will at first continue to increase with the increasing volume 
intercepted by the spheres, but will not do so right up to 
5 + i/zp, because the intensity in the rod is greatest in 
the middle and falls off at the outside. The maximum on 
the film will therefore lie somewhere between the value of 
corresponding to S = S and that corresponding to 
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S = S + 1/2/), or at about S + i/^p, say. The peak of 
tlTe reflection is therefore displaced from the position where 
it would be expected by 

AO - (A/2 cos 0)( 1 /4/>) = o- 1 25 X/p cos 0. . (3 5) 

This displacement was discovered experimentally before 
the peculiar nature of the reflections was realized, and was 
interpreted as a contraction of the lattice in small crystals 
of graphite. 

The approximate line profile as a function of s is shown 
in figure 6. The concept of integral breadth is not useful 
for dealing with such reflections, as the intensity never 
falls quite to zero with increasing s, and it is necessary to 







"Vp O '/p Vp s 

FlG. 6. — Approximate line profile for hkl reflections. 

use the half breadth as a measure of the broadness of the 
reflection. The intensity reaches its maximum for 
s ~ i/^p, but from the volumes intercepted by the 
spheres in the two positions it is pretty clear that the inten- 
sity for s = o is greater than half that for s = 1/4/). At 
a guess we might say that the intensity is half the maximum 
at the same distance to the other side of 5 = o, or at 
s = - i/4/». The distance on the low-angle side between 
the maximum and the half value is thus i/2p. On the 
high-angle side there is nothing to aid in the estimation 
of the distance it is necessary to go to reach the half value, 
but as the rise in intensity is steeper than the fall it will 
certainly be more than on the low-angle side, and as a guess 
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we might double it and say i/p. The total half breadth in 
s of the reflection would then be 3/2^, and the half breadth 
in 20 would be 

ft' = (A/cos 0)( 3 / 2 p) = i-sl/p cos 0, . (36) 

and the apparent particle size calculated from the half 
breadth would be 



e' = X/P' cos = o-6yp. 



(37) 



If our estimates are anywhere near right, therefore, the 
numerical factor is considerably different from the 1-13 
derived above for the 00/ reflections. Warren (1942) has 
calculated more exact values, assuming that the intensity 
falls off exponentially with the square of the distance from 
the centre of the rod, and obtaining 



and 



£ ' = o-57p . 
AO =o-i6X/pcos6, 



(38) 
(39) 



so that the rough estimates made here give a semi-quanti- 
tative explanation of the broadening and displacement of 
the lines. The apparent particle size can be estimated 
from both the apparent particle size and the displacement. 
4. Line broadening due to three-dimensional mistakes 
in AuCu 3 has been extensively investigated by Jones and 
Sykes (Sykes and Jones, 1936 ; Jones and Sykcs, 1938), 
who have measured the apparent particle size for various 
reflections. These are given for a typical specimen in 
table II, most of the entries being average values from four 
or five photographs. It will be seen that e varies markedly 
with h, k and /, the extreme values, those for 100 and no, 
differing by a factor of two. The reason for the broaden- 
ing is fairly clearly the fact that the type of ordering is not 
the same throughout each single crystal, but ' changes 
step ' occasionally, so that the probability of two cells 
separated by a distance t in the h k I direction being the 
same decreases with increasing t. The boundary between 
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regions of different types of order must be plane, as it 
would othenvise be unstable (Bragg, 1940), but there is no 
other general restriction on the mode of occurrence of 
mistakes. Wilson (1943) has considered five different 
possibilities : 

(i) mistakes occurring on any plane, whatever its 

inclination to the crystal axes, 
(ii) mistakes occurring only on planes perpendicular to 

the edges of the unit cube «ioo> planes), 
(iii) mistakes occurring only on planes perpendicular to 

a face diagonal of the unit cube «no> planes), 
(iv) mistakes occurring only on planes perpendicular to 
the body diagonal of the unit cube «i 1 1> planes), 
and 
(v) mistakes occurring in such a way that gold atoms do 
not come into contact. 
The fifth of these resembles the second in some ways, as 
a little consideration shows that it is only on planes per- 
pendicular to the cube edges that mistakes can occur with- 

TABLE II 
Observed and calculated apparent particle sizes in a specimen 
Audi, 



Line 


fob,. 


tub. 


% error 


100 


200A 


203A 


+ i-5 


no 


IOO 


109 


4- 9-0 


210 


120 


Hi 


+ 9-2 


211 


142 


135 


-5'2 


221 


"5 


112 


-2-6 


310 


122 


114 


-6-6 


320 


150 


152 


+ 1-2 


321 


124 


122 


- 16 


410 + 322 


128 


127 


-o-8 


411 


161 


>57 


- 2-5 


421 


"5 


1 12 


-26 



out bringing gold atoms into contact, but there is nn added 
restriction. If, for example, the order is of type A, and 
a mistake occurs on a plane perpendicular to the a axis» 
on the second hypothesis the order could change to B of 
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C or D, while on the fifth it could change only to D. The 
complete list of mistakes possible without bringing gold 
atoms into contact is 

on planes perpendicular to the a axis, A ^ D, B ^ C, 
on planes perpendicular to the b axis, B ^ D, A ^ C, and 
on planes perpendicular to the c axis, C ^ D, A * B. 
Values of J(t) for each of these modes of occurrence of 
mistakes can be worked out by methods similar to those 
already employed for layer structures, but slightly more 
complex, because of the three-dimensional nature of the 
mistakes. If h, k and / are all even or all odd 

7(0 = (Au + 3 Cu)*, 
where Au and Cu are the scattering factors of the gold and 
copper atoms respectively. These are the reflections 
common to the ordered and the disordered state, and are 
sharp, in agreement with experiment, since J(t) is constant. 
If h y k and / are not all even or all odd the five possibilities 
of occurrence of mistakes lead respectively to the following 
values of /(/) and e : 

Variety of Apparent 

Mistake J{t) Particle Size 

(i) (Au- Cu) s exp{-40tU|/3a) 30/20. 

(ii) (Au- Cu) a exp<- 4 «A|<l/3 VJVa} 3VNa/zal, 

(iii) (Au-Cu)-exp{- 4 a(/i+A)|f|/ 3 v' (zAf)o) 3 V(2A0fl/2a(>/+ k) 

(iv) (Au-Cu)-'exp{-4a(/i+/ ! +/)|/|/3V , (3N)fl} 3 V(3N)a/2 a(A+A+/) 

(v) (Au - Cu) 2 cxp{- 2 a(/> + q)\t\/VNa } y/Na/tx(j> + q), 

where h, k, I are the positive values of the indices arranged 
in the order h > k > /, N = h 2 + k- + I 2 , a is the chance 
of a mistake occurring in a distance equal to a (the edge 
of the unit cell), and p and q are the indices of the same 
parity (that is, the two that are even, or the two that are 
odd, as the case may be). It is not possible to choose a 
so that any of these expressions gives quantitative agree- 
ment with the experimental values reproduced in table II, 
but the fifth comes close to it. The chief discrepancy is 
the 100 reflection, for which the theoretical value of e is 
Q /a(o + o) = co, so that the reflection should be sharp. 
It is actually the sharpest of the broadened reflections, 
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which suggests that a combination of two varieties of 
mistake, the fifth which gives good qualitative agreement, 
and the first to give a general fuzziness independent of 
indices, would give quantitative agreement. For this 
combination it is found that 



(vi) /(/) = (Au - Cu) 2 exp { - [ 2 a'(/> + q)/y/N 

+ 4a73](M/«)}. • • • (40) 

where a' and a" are the probabilities of the two kinds of 
This equation gives 

• • (4i) 



mistake. 



e = 



a'(/> + q) + 2a; 



■\/Na 2 a 

With this expression for e it is possible to get fair 
quantitative agreement with the experimental values. 
Apparent particle sizes calculated with a.' /a = 0-00303 
and 2a" /3a = 0-00492 are given in column three of 
table II, and it will be seen that the greatest difference 
between the observed and calculated values is less than 
10 per cent. This agreement is very satisfactory, the 
difficulties of measuring line breadths being considered, 
but it does not seem impossible that some other mode of 
occurrence of mistakes would give even better agreement. 
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CHAPTER VII 

THE USE OF COMPLEX QUANTITIES 

1. In the earlier part of this book an attempt has been 
made to keep the necessary mathematics as simple as pos- 
sible, but a stage has now been reached where it seems 
desirable to introduce complex quantities. It is possible 
to proceed a little further without making explicit use of 
them, but already in the discussion of crystals with mis- 
takes it has been necessary to postpone one or two topics 
because elementary methods became very clumsy, and 
clumsiness is evident from the first in any attempt to dis- 
cuss diffraction by distorted crystals. It is, of course, 
impossible to give an adequate discussion of the nature 
and properties of complex quantities in the space available, 
but it is convenient to summarize the properties that will 
be used in the following chapters, and incidentally to clear 
up the difficulties postponed from chapters V and VI. 
A complex quantity is of the form 

z = x + iy, . . . (1) 

where i- = - 1 and x and y are ordinary numbers or vari- 
ables. It is usual to call x the real part of z and y the 
imaginary part. A complex quantity may be represented 
in a plane by a line drawn from the origin of length equal 
to r, the positive value of V(* 2 + y 2 )> and making an angle 
with the horizontal axis, where r cos = x, r sin =y. 
Since 

exp(«0) = 1 + i0 -0 2 /2! -z'0 3 /3! +074! +/O75! - — 

= 1 - 072! + 074! - — + i(0 - 073! + 075! - •••) 
= cos + i sin 0, . . . . . (2) 

a complex quantity may be written 






z = x + iy 
= r cos + ir sin 
= r exp (iO). . 
85 



(3) 



86 



X-RAY OPTICS 



. 



The sum of two complex quantities z y and z., is 

ft, + z 2 = Xl + iy x +x.,+ iy 2 m Xl + X% + i( yi + y 2 ), . ( 4 j 

and their product is 

*i*2 = ( x x + iyi)(*t + iy z ) 

= (%*, -y 1 y 2 ) + t'(^y 2 + xsyj. 

The product may be written more neatly in terms of r 
and 0, for 

z x z z = r, exp (t'0,) .r 2 exp (i0 a ) 

= r 1 r 2 expO(0 1 +0 2 )}. . . (6) 

If the sign of i is changed wherever it occurs in a complex 
expression, the result is called the conjugate of the original 
expression. It is usually represented by the same letter 
with an asterisk affixed, so that if z = x + iy, 

z* = x -iy. . . . ( 7 ) 

A quantity multiplied by its conjugate is purely real, and 
equal to the square of its absolute magnitude, for 



zz* = (x + iy)(x - iy) = x 
The relations between z, x, y, r, 



+ y* = r* 



• (8) 
and z* are shown 
diagrammatically in 
figure i for two values 
of z. In forming the 
conjugate of a com- 
plex expression it is 
not necessary to re- 
duce it to the form 
x + iy first. For in- 
stance, the conjugate 
of (A - iB)/(C + iD) 
is (A + iB)/(C - iD)- 
2. These properties 
of complex quantities make possible a compact expression 
for the amplitude scattered by a unit cell. This is repre- 
sented both in amplitude and phase by 

F =Zfj exp {zmihxj + hy i + ty)} . (9) 



z 


» 


Y 




7 






x (] 


><0 x 




•r 




V s 


/r 


c$W 






X 


z 






/ 


2 


► 



Fig. i. — Representation of complex 
quantities. 
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where fj is the scattering factor and g* y^ Zj arc the co- 
ordinates of the y lh atom in the unit cell, and •■ means that 

the expression is to be summed for all atoms in the unit 
cell. It is easily verified that the real part of this is equal 
to A and the imaginary part to B (chapter V, paragraph 1). 

2. In the calculation of the scattering factor from various 
imperfect structures it is necessary to calculate the mean 
value J of AA' + BB' and K of AB' - A'B, where A, B 
and A', B' are the A and B parts of the structure ampli- 
tudes of two unit cells with a certain separation in the 
crystal. If F = A + iB, F' = A' + iB', 

FF'* = (A + iB){A' - iB') 

= (AA' + BB') - i(AB' - A'B). . (10) 

Thus J is the mean value of the real part of FF'* and - K 
is the mean value of the imaginary part. It is often con- 
venient to represent / - iK by the single letter Y. 

3. It is therefore important to calculate the mean value 
of expressions like exp (27tiux), where u is fixed and x is 
some variable quantity, such as one of the co-ordinates of 
position within a crystal, or the displacement of a unit cell 
from its mean position. When x varies over a small range 
of values on either side of zero, an approximate mean value 
can be obtained by expanding the exponential in a power 
scries : 

exp (zTiiux) m 1 + zniux - znWx 1 - piHuW/i + •••. 

If positive and negative values of x are equally frequent, 
the mean values of x, x 3 and other odd powers of x will 
be zero, and 

< exp (zniux)'} = 1 - 2te 2 w 2 <* 2 > + 2ti; 4 h 4 <# 4 >/3 - •••, 

where the angle brackets indicate average values. Now 

exp (- 2?r 2 tt 2 2 » = 1 - 27ihi\x*) + 2?r 4 M 4 <* 2 > 2 - -, 

s o that, if the difference between <* 4 >/3 and <x 2 > 2 can 
°e ignored, 

< exp (zniux)} = exp (- 27t 2 M 2 <^ 2 ». 
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In certain cases this approximation is exact, but it is always 
fairly good for small values of x. 

4. If the probability of x having any particular value is 
given byf(x), where the integral oif(x) from - 00 to + 00 
is unity, the mean value of cxp (27tiux) is given by thi 
integral 




<exp (27iiux)y 



J 



f(x) cxp {zniux) dx, 






which will of course depend on u. For the present, how- 
ever, let us forget about the interpretation of this expression 
as an average, and consider the function g(u) related to 

fix) by 

g(u) = $f(x) exp (zniux) dx, . . (n) 

where the limits are understood to be - 00 to +00 and 
f(x) is any function of x for which the integral exists, and 
does not necessarily fulfil the condition that jf(x) dx = 1. 
The function g(u) so defined is called the Fourier transform 
of /(*). It has the curious property that 

f( x ) ~ fg( u ) ex P (~ 2nixu) du. . . (12) 
Except for a change in the sign of i,f(x) is also the Fourier 
transform of g(u). If both / and g are purely real they are 
Fourier transforms of one another, without this slight lack 
of symmetry in the relation. There have been numerous 
examples of somewhat-disguised Fourier transforms in 
earlier chapters ; for example, the amplitude of reflection 
from a small crystal (chapter III, paragraph 6) is the three- 
dimensional Fourier transform of a function with value 
unity inside the crystal and value zero outside, and the 
intensity of reflection from a layer structure with mistakes 
(chapter V, equation (29) ) is the transform of /,„ - iK m . 

5. A rigorous proof of the relation expressed by equa- 
tion (12) is difficult, but depends on the following argument 
From the definition 

g{u)=jf(t)exp{2mut)dt, . . (13) 
since the substitution of t for x as the variable of integration 
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makes no difference. Let 

F(x) = fg(u) exp (- zm'xu) du. . . . (14) 
Then 

F(x) = JLT/(0 ex P (27tiut) dt] exp ( - 27tixu) du, 
or, on changing the order of integration, 

F(x) - /Lf ex P M< " *)"} du]f(t) dt. . . (15) 
The inner integral, J exp {zjii{t - x)u} du, is very peculiar. 
If t = x it has the value infinity, since exp (o) is unity, 
whereas if t^x exp {zm{t - x)u) continually oscillates, but 
has the average value zero. The inner integral is thus zero 
if t^x, but at t = x it has an infinitely high and infinitely 
narrow peak. In the integration with respect to t, there- 
fore, the only region contributing anything will be an infi- 
nitely narrow strip containing the point t = x, and F(x) 
will be equal tof(x) multiplied by the area under the peak. 
In various ways this can be shown to be unity, so that 

F(x)=f(x) . . (16) 

The difficulties over which much trouble has to be taken 
are first, justifying the change in the order of integration, 
and secondly, showing that the area of the peak is actually 
unity. One method of approach is to consider 

F(x, X) = fg(u) exp ( - hi* - znixu) du . (17) 

instead of F(x), and to show that F(x, A) tends to f(x) as 
X tends to zero. 

6. Let us now return to the problem of the scattering 
of X-rays from an imperfect crystal. The rays scattered 
by the unit cell at the point n x a, nj), n 3 c (the n th unit cell, 
for short) are advanced in phase with respect to those 
scattered by the cell at the origin by 

<5 n = zn^ftjtt + n 2 v + 7i z w) . . (18) 

Where u, v, w are the co-ordinates of Q with respect to one 
of the points of the reciprocal lattice (chapter III, para- 
graphs 1, 6). If F n is the value of the structure amplitude 
(9) for the » th unit cell, the amplitude scattered by the 

7 
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whole crystal will be 

G(u, v, to) = EF n exp («5 n ), 



■ (19) 



where Z means that the expression is to be summed for 

n 

all unit cells. The intensity scattered will be given by 
equation (19) multiplied by its complex conjugate, so that 

H(u, v, w) = ZZF n F;, exp {/(<5,, - d,,)} 

n if 

= £ZF n F* n > exp {znilu^ - n\) + v{n, - n'.,) 

+ w(« 3 -n 3 )]}. . (20) 

Note that in calculating H it is necessary to sum for all 
unit cells separately in G and G* ; 

/ = £F n F: exp {i(d n - (5,,)} = r|F n |* . ( 2I ) 

n ,„ 

is something quite different from H, since it is independent 
of u, v, w. It is in fact the total intensity scattered by the 
crystal in the hkl direction, or the integral of // over one 
unit cell of reciprocal space. 

In order to obtain from (20) an equation suitable for the 
actual calculation of H, it is necessary to get rid of the 
separate values of F for each unit cell, winch are not known 
individually, and replace them by some sort of average 
values, which can be calculated in ways similar to those 
used for layer structures. Clearly the average must be 
taken in such a way that the exponential factor in equa- 
tion (20) is not changed, or in other words, F n F' n > must be 
averaged for all values of n and ri that give the same values 
of 0J - ti\, n 2 - 11'.,, « 3 - n' a . In equation (20) let 

n\ m m + tn lt n'., = » 2 + m 2 , ri a = n 3 + m 3 , . (22) 
so that 

H(u, v, zo) =ZZF n F' l 4„ l exp { -zni^u +m 2 v +m a w) } 

n m 

=ZIZF n Fn +m \ exp { -zxiimjii +m 2 v +m. 6 w) }, (23) 
proper attention being paid to the limits of the summations- 
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What does the sum in brackets mean ? Starting with any 
cell, designated by the numbers n v n 2 , n a , we are to find the 
cell m x a away in the direction of the a axis, m 2 b away in the 
direction of the b axis, and m 3 c away in the direction of 
the c axis, then multiply the structure amplitude of the 
first cell by the conjugate of the structure amplitude of 
the second, and add for all possible unit cells. For how 
many cells is it possible ? At first one might answer : for 
every cell in the crystal, but actually this is true only for 
m i " m 2 = m z = °- F° r % — 1, for example, a layer of 
cells all over one side of the crystal must be omitted, as 
there are no cells at a distance a away in the direction of 
the a axis. For m x = io, a layer ten cells thick must be 
omitted. Consideration shows that in general the number 
of cells to be included in the summation is the number 
common to the crystal and its ' ghost ' shifted a distance 
m r a in the direction of the a axis, m 2 b in the direction of 
the b axis, and w 3 c in the direction of the c axis. If this 
volume is called V m , the number of cells to be included 
is V m /U, where U is the volume of one unit cell. If, 
further, Y m is the average value of the product F n F n+m 
for all unit cells related by the displacement w x a, mj), m 3 c, 
then equation (23) becomes 

H{u, v, w) = U- l ZV m Y m exp { -2jri(;«i« +m 2 v +m s w)}. (24) 
m 

This is the general expression for the intensity of reflection 

as a function of position in reciprocal space for a crystal 

whose unit cells arc not all alike, and is the solution of the 

problem postponed in chapter V, paragraph 9. It is clear 

that V m is a generalization of the V(t) used in chapter IV, 

and Y m of J m - iK m , where J m and K m are the functions 

defined in chapter V. From the definitions it follows that 

changing the signs of the m's leaves V. m unchanged, but 

turns Y m into its conjugate. A little care is necessary, 

therefore, in carrying out calculations, as the algebraic 

expressions for V m and Y m do not always take account of 

these properties automatically. 
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The summation over m in equation (24) is to be done 
for all values of m^, m 2 , m 3 for which the product to be 
summed exists. Since the sign of the non-existence of the 
product is the vanishing of V m , \vc may write the limits 
symbolically as ± co for each of the m's, remembering to 
stop the summation when V m vanishes. If, u, v, w are 
not too big, and Y m is a ' smooth ' function of the m's, 
the summations in equation (24) may be replaced by 
integrations, giving 

/•+» r+x r+K 

H{u,v,w)={Uabc)-A V(x,y,z)Y(x,y, 

J —to J —to J — X 

exp {- 27ii(ux/a + vy/b + wz/c)} dx dy dz, . (25) 

where x = m x a, y = m 2 b, z = m 3 c, and V(x, y, z), Y(x, y, z) 
are V„„ Y m with the appropriate substitutions. In a crystal 
containing mistakes Y(x, y, z) may decrease so rapidly that 
the variation of V(x, y, z) can be ignored, and the intensity 
written 

H(u, v, w) - (Uabc)-WJjjY(x,y t z) 

exp {- 27ii(ux/a + vy/b + wz/c)} dx dy dz. . (26) 

For a crystal without mistakes Y(x, y, z) is a constan 
equal to F 2 , and the intensity becomes 
H(u, v, to) = (Uabc)- l F°iHV(x, y, z) 

exp {- 27ti(ux/a + vy/b + wz/c)} dx dy dz, . (2 
giving an expression alternative to equations (15H17) of 
chapter III. It is easy to verify that equation (27) gives 
the same results in particular cases as were obtained by 
elementary means. For instance, for a parallelepiped of 
edges X, Y, Z, 

V{x, y, z) = (U/abc)(X - \x\ )( Y - |;y| )(Z - |*| ), (28) 
and equation (27) gives on integration 
H(u, v, w) 

(29) 









p2 sin 2 (nuX/a) sin 2 (nvY/b) sin 2 (nwZ/c) 



(nuY {nvY ' (nwf 

in agreement with chapter III, equation (7). 
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7. Equation (26) may be used to discuss the intensity of 
reflection from a structure with a variable interlayer spacing, 
postponed from chapter VI. The values of /„, and K m 
found there (end of paragraph 1) give 

*« = Jm ~ ***-m 

= F 2 [(i - a)" 1 + »ia(i - a)'"- 1 exp (- znilg) 

+ m <P ~ I ) a 2( I _ a )».- 2 exp ( _ 27li . 2 ig) 

21 

+ ••• + a m exp (- 27iimlg)] 
= F 2 [i - a + a exp (- 27tilg)] m . . . (30) 

by the binomial theorem. The m in this equation is what 
would be called m 3 in the notation of the present chapter ; 
since Y m depends only on w 3 and not on m 1 and m 2 , 
H{u, v, w) drops rapidly to zero very rapidly with u and 
w. If a is small 1 - <x[i - exp (- znilg)] is approximately 
equal to exp {- a[i - exp (- 27iilg)]}, so that 

Y m '=. F 2 exp {- w 3 a[i - exp (- 27zilg)]}, 
Y(z) = F 2 exp {- a[i - exp (- 27iilg)]z/c}. . (31) 

Equation (26) for the intensity of reflection as a function 
of u, v, w reduces to 

H{p) = {V/U) [ + * Y(z) exp (- 27tiwz/c)d(z/c), (32) 

J —00 

or, since Y(- z) = Y*(z) (page 91) and V/U = N, 
H(w) 

= NF 2 |f exp { [ - 27tiw -a { 1 - exp ( -zmlg) }]z/c }d(z/c) 

+ complex conjugate 

= NF z [{ol[i- exp (~27iilg)] + 2711'w}- 1 

+ complex conjugate] 
= JVF 2 [{a(i - cos 2nlg) + i{27iw+v. sin ZTilg)}- 1 

+ complex conjugate] 

=iV F 2 2a(i-cos2^ ) 

a 2 (i - cos 27ri^) 2 + {zrvw+u. sin 27r^) 2 ' 
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Thus the maximum intensity of reflection is displaced by 
(a/27t) sin znlg towards the origin of reciprocal space, and 
the integral breadth in to is (a/2)(i - cos znlg). The dis- 
placement is in the direction to be expected, for the mistake 
increases the average c spacing, and hence decreases the 
distance of the maximum intensity from the origin of 
reciprocal space. Both the displacement and the broaden- 
ing depend on /, and for reflections such that Ig is integral 
they vanish, and become very small if Ig is nearly integral. 
The X-ray reflections will therefore be affected to very 
different extents, depending on the value of /. Hendricks 
and Teller (1942) have calculated in detail the 00/ reflec- 
tions of a partially dehydrated montmorillonite, which 
appears to consist of layers of thickness 10 A and 13 A 
with about equal frequency and a about one-half. They 
give no quantitative comparison with experiment, but state 
that the qualitative agreement is good. Equation (33) 
cannot be expected to hold for such a large value of a, 
but, since Y m is of the form p m , H(w) can be got by direct 
summation of equation (24) as a geometric series. 

8. The methods of chapter V, paragraph 9, were suffi- 
cient for expressing the profiles and integral breadths of 
lines on powder photographs of imperfect crystals. Never- 
theless, expressions in terms of complex quantities are some- 
times more convenient for calculation. If Y(t) =J(t) -iK(t) 
is the mean value of the product FF* of the structure 
factors of two cells separated by a distance t in the hkl 
direction, it is readily found that the line profile is given by 
I{S) = U-^V{t)Y{t) exp (27tist) dt, . . (34) 
the apparent particle size is, if the maximum is not displaced, 
e=I s /I = [VY(o)]-*IV(t)Y(t)dt, . (35) 
and the integral breadth is given by X/e cos 0. For 
crystals in which Y(t) decreases rapidly with t, so that the 
variation of V(t) can be neglected, these reduce to 

I(S) = (V/U)jY(t) exp (2711st) dt . . (36) 
and e-[Y{o)]-*fY(t)dt. . . . (37) 
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9. Equation (26) may be used to discuss the distribution 
of intensity in reciprocal space for single crystals of AuCu 3 
containing mistakes. This was of course unnecessary in treat- 
ing line broadening in Jones and Sykes' powder photographs, 
but the actual distribution in reciprocal space can be investi- 
gated by the scattering from a single crystal (Guinier, 1945 ; 
MacGillavry and Strijk, 1946). Sufficient details have not 
yet been published to show whether or not the distribution 
of mistakes postulated to explain the powder photographs 
is correct, but it is of some interest to obtain expressions 
for the intensity of scattering for the various models. 

The values of Y{x,y, z) required follow readily from the 
values of J(t) for AuCu 3 given on page 83. Since /(f) 16 
the mean value of AA' + BB' for cells separated a distance 
/ in the hkl direction, the real part of Y(x, y, z) will be /(/) 
with ht/\/N replaced by x, kt/y/N by y and lt/y/N by z. 
The imaginary part is zero, since K{t) is zero. Thus for 
the five varieties of mistake Y(x, y, z) has the values 

(i) (Au - Cu) 2 exp { - 4 aV(* 2 + y 2 + ^)/ 3 a}, 

(ii) (Au-Cu) 2 exp{- 4 a|*'|/3«}. 

(iii) (Au - Cu)* exp {" 4*( M + l/l )/3Vm). 

(iv) (Au - Cu) 2 exp {- 4 a( |x| + \y\ + \z\ )/3V3*}> 

(v) (Au - Cu) 2 exp {- 2 a( |*"| + \y"\ )/«}. 

In expressions (ii) and (iii) x', y', z' represent x, y, z 
arranged in the order \x'\ > \y'\ > \z'\, and in expression (v) 
x" andjy" are the co-ordinates corresponding to the indices 
of the same parity. The combination of mistakes of the 
first and fifth kinds gives 

(vi) (Au - Cu) 2 exp {-[2a'( |*"| + \y"\ ) 

+ 4 a'V(* 2 +>' 2 +^)/3]/«}. 
where a' and a" are the respective probabilities of the two 
kinds of mistake. For a cubic crystal equation (26) becomes 

/•+» /.+» -.+«> 
H(u, v, to) = Na-*\ Y(x,y,z) 

J —x J — « J — » 

exp { - 27ii(ux + vy + toz)/a } dx dy dz, (38) 
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and the problem is now only that of performing the inte- 
grations with the values of Y(x,y, z) given by expressions 
(i)-(vi). This is perhaps easiest for (iv) : 

H(u, v, w) 

- iV(Au-Cu)^j +X exp {-[ 4 a \x\ / Zy / 3 

+ 2mux]/a}d(x/a) 
multiplied by similar integrals in y and z 

= iVtAu-Cu^fTexp {-[4x/ 3 V3+zrc™]{x/")}d(x/a) 

+ complex conjugate 

multiplied by similar expressions in y and z 

= A^Au-Cu) 2 !" "__i r + ' _1 

_4a/3 V3 +2mu 41X./3 V3 -2a»«J 

multiplied by similar expressions in y and z 

= AT(Au-Cu) 2 

8a 8a 8a 

3V3 3V3 3V3 



teh^* <&+» (wih™r 

The integration for x" and y" in (v) is similar, giving 
H(u, v, w) 



(39) 



= N(Au - Cu) 2 - 



jjr,.F(t0") (40) 



a 2 + (ttu") 2 a 2 + {tw") z 

where u" and v" are the variables in reciprocal space cor- 
responding to the indices of the same parity, and F(w") 
is the ' normal ' variation of intensity of reflection with 
the third variable ; it is large only for w" ~ o, and falls 
to inappreciable values for w" of the order of the lattice 
parameter divided by the particle dimensions. This is a 
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very peculiar distribution of intensity in reciprocal space ; 
the regions of high intensity are disc-like, with the plane 
of the disc parallel to the directions indicated by the indices 
of the same parity. Thus for the 100 reflections the disc 
is perpendicular to the 100 direction, for the 010 reflection 
to the 010 direction, for the no reflection to the 001 direc- 
tion, for the 211 reflection to the 100 direction, and so 
forth. The distribution of intensity in reciprocal space is 




Fig. 2. — Reciprocal lattice for AuCu 3 . 

represented diagrammatically in figure 2. The calculation 
for expression (vi) is more difficult, but it seems that the 
result would be qualitatively similar to that with (v). The 
regions of high intensity would be discs with the same 
orientations, but of appreciable thickness. Unfortunately 
the single photograph reproduced by Guinier (1945) is 
insufficient to show if the regions of high intensity are of 
this form or not. It does show, however, that the sharp 
reflections are surrounded by a region of diffuse scattering 
that cannot be explained by mistakes of the type treated 
here. It is probably an effect of lattice distortion, the 
copper atoms surrounding a gold atom being displaced 
somewhat from the lattice points. 

The integration for expression (i) is most easily done by 
transforming to polar co-ordinates with the axis passing 
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through the point u, v, w. If s = \/(u 2 + v 2 + w 2 ) and 
r = \/(x 2 + y 2 + z 2 ), equation (38) becomes 

H(t) =N\ f f exp {-[4a/3 + zjcis cos 9](r/ a)} 
Jo J J 

(r/a) 2 sin Od(r/a) dd d<j> 
..f" 1 2sin0<f0 

= 27lN 






[4a/3 + ZTEWcosfl] 3 
8a 



= 4?cN r 



A3 



"4- <*»)«]' 



(4i) 



As would be expected, the distribution is spherically sym- 
metrical, falling off ultimately as N/qnW from its initial 
maximum of 2j7tN/8ix. 3 . The expressions for H corre- 
sponding to (ii) and (iii) can be obtained similarly, but 
are of considerable complexity, and the simpler cases 
illustrate the methods sufficiently well. An experimental 
investigation of the scattering in the neighbourhood of 
several points of the reciprocal lattice to show definitely 
if any of the models is correct would be of great interest. 
A preliminary investigation by Edmunds, Hinde and Lip- 
son (1947) and the work of Strijk and MacGillavry (1946) 
on reflections of the type hkli are in qualitative agreement 
with model (vi). 
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CHAPTER VIII 
DISTORTED CRYSTALS 

1. In the crystals discussed so far the lattice has been 
perfect, or, in other words, the unit cells have been all of 
the same size and shape, though the contents of the cells 
may have varied throughout the crystal. Crystals in which 
the size and shape of the cell vary from place to place form 
an equally important but more difficult problem. The 
distortion of the lattice may be of several kinds. For 
example, at all temperatures above the absolute zero the 
atoms of a crystal are kept in continual thermal agitation. 
The agitation is best imagined as the superposition of 
a number of elastic waves traversing the crystal. These 
waves are similar to sound waves, but the amplitude is in 
general smaller and the frequency higher. The frequency 
is, however, low compared with that of the X-rays, and 
for purposes of calculation the atoms may be considered 
at rest, but displaced from their mean positions. This 
problem is treated in the next chapter. Other sorts of 
distortion are found in metals used for constructional pur- 
poses. Perfect crystals of most metals are soft and easily 
deformed, and metals in a useful form are generally in 
a state of internal strain. This may be produced mechani- 
cally in drawing or rolling, or ' chemically ', by treatment 
of an alloy in such a way that it begins to separate into two 
phases, but is not permitted to finish the process. Thus 
duralumin (about 96 per cent aluminium, 4 per cent 
copper) in equilibrium at room temperature would consist 
of a mixture of crystals of practically pure aluminium and 
of a compound Al 2 Cu. As used, however, Cu or Al 2 Cu 
separates in sheets within the aluminium crystals, but on 
lattices continuous with those of the parent crystals. The 
unit cells of the aluminium and the new structure will not 
fit together without a considerable amount of strain, and 

99 



IOO 



X-RAY OPTICS 



the composite crystal is much harder than pure aluminium. 
In a sense the Cu or Al a Cu sheets may be regarded as 
mistakes, but there is a variation of lattice parameter as 
well as of cell contents. Another example of a variation 
of lattice parameter within the individual crystals is given 
by alloys with compositions in the neighbourhood of 
Cu 4 FeNi 3 . The alloys are in themselves of no particular 
importance, but it happens that the lattice parameter is 
approximately a sinusoidal function of position within the 
crystal, so that the theoretical treatment is comparatively 
straightforward and the interpretation of the photographs 
easy. 

2. The alloy Cu 4 FeNi 3 has been investigated by Daniel 
and Lipson (1943, 1944). In equilibrium above about 
800 c ' C. it consists of cubic crystals of one kind only, but 
below this temperature it consists of a mixture of cubic 
crystals of two kinds, differing in composition. If a speci- 
men in equilibrium above 8oo° C. is held at a lower tem- 
perature (say 650 C.) for a time insufficient to reach 
equilibrium, a state intermediate between that of perfect 
crystals of one kind and that of perfect crystals of two kinds 
is obtained, and can be maintained indefinitely by cooling 
the specimen to room temperature. In the intermediate 
state the lines on powder photographs are sharp, but are 
flanked by satellite lines, whose existence can be explained 
by the presence in the crystals of layers perpendicular to 
the crystal axes and alternately copper-rich and copper- 
poor. These will have a varying lattice parameter in the 
direction perpendicular to the layers, and the assumption 
that the variation is sinusoidal accounts for most of the 
details of the photographs. Since Fe, Ni and Cu have 
approximately the same atomic number the structure ampli- 
tudes of the cells will vary little, and may be taken as 
independent of position. 

Consider, therefore, a crystal approximately cubic, but 
with the lattice parameter varying sinusoidally in the 
x direction from a - S to a + <5, the wave-length of the 
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variation being L. Then the position of the n th unit cell 
will be given by 



X n = n x a + (Ld/zna) cos (znn^/L), 



(0 



since the distance in the x direction between the origins of 
successive cells is 

[(«! + i)a + (Ld/zna) cos {zn(n x + i)a/L}] 

- [n^a + (Ld/zna) cos (znn^/L)] 

= a + (Z^/27ra)[cos (znn^a/L) cos (zna/L) 

- sin (znn x a/L) sin (zna/L) - cos (znn^/L)] 

which becomes 

a - 6 sin (znn^a/L), . . (2) 

since a is small compared with L, and cos (zna/L) "— 1, 
sin (zna/L) '—. zna/L. If h', k', V (not necessarily integral) 
are co-ordinates in reciprocal space, the amplitude of 
reflection as a function of position in reciprocal space will 
be given by 

G(h', k', V) = FZexp {zni[h'x H + k'y n + l'z n ]/a} 

n 

= FE exp \zni\nji' + (Ldh' /zna 2 ) cos (znn x a/L) 

+ n 2 k'+n z l']} (3) 

where F is the structure amplitude of one unit cell. As 
a function of k' and /' this is normal, and has large values 
only for k' = k, I' = /, where k and / are integral. Its 
behaviour as a function of h', however, requires special 
investigation. Let us consider therefore 

G(h') = FE cxp {zni[nh' + (Ldh' /zna 2 ) cos (znna/L)] } 

= FEexp {znink'}.exp {zni(L6h' /zna 2 ) cos (znna/L)} (4) 

n 

If 6 is so small that Ldh' /a 2 is small compared with unity, 
the second exponential in equation (4) may be replaced by 
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I0 3 



the first two terms of its power series : 
exp {27ti(L61i /271a 2 ) cos {zTcna/L)} 

= 1 + i(L6h'/a 2 ) cos (zTtna/L) 
- 1 + i{L6h' /2a 2 ) exp {zjcina/L} 

+ i(Ldh'/2a 2 ) exp {- 27iina/L}, . (5) 
and equation (4) becomes 
G{h') = FZexp {27tinh'} 

" + iF{Ldh' /2a 2 )Z zxy {zninW + a/L)} 

+ iF(Ldh'/2a 2 )Zcxp {ito(// - a/L)}. . (6) 
n 

The first of these sums has a maximum of amplitude pro- 
portional to F at h' = h, where h is any integer ; the second 
has a maximum of amplitude proportional to F{L6h' /2a 2 ) 
at h' + a/L = h, where h is any integer ; and the third has 
a similar maximum at h' - a/L = h, where h is any integer. 
Thus, in addition to a strong reflection at the points of the 
reciprocal lattice, there are weaker ones at distances of 
± a/L from it. The intensities are in the ratio 

[Ld(h - a/L)/ 2 a 2 ] 2 : 1 : [L6(h + a/L)/za 2 ] 2 , 
or approximately (IJih/za 2 ) 2 : 1 : {Ldh/za 2 ) 2 . The inten- 
sities of the extra maxima, or ' side-bands ' as they have 
been called, increase as the square of the order of reflection. 
3. For comparison with this result it is interesting to 
work out the reflection from a structure in which the 
structure amplitude, instead of the lattice parameter, varies 
sinusoidally. 

If *»-»«. \ tn\ 

F n = F[i + gcos{27tna/L)l] ' - w 

where g is a fraction and L is the wave-length of the 
variation, the amplitude of reflection is 
G(h') =FZ[i +g cos (27tna/L)] exp {zninh'} 

-.FT exp {27tinh'}+F(g/2)Zexp temdk'+a/L)\ 

n n 

+F(g/2)Zexp tenMh'-a/L)}. (8) 

n 

As before, the three sums lead to maxima at h' ~ h, 



h ± a/L, where h is any integer, but the intensities are in 
the ratio g 2 /$ : 1 : g 2 /^, so that the side-bands are the same 
in all orders. In particular, there are side-bands flanking 
the origin of the reciprocal lattice (h = o), whereas these 
are of vanishingly small intensity for lattice parameter 
variation. Bragg and Lipson (1943) have devised gratings 
that illustrate the different diffraction phenomena caused 
by mistakes, variation of structure amplitude, and variation 
of lattice parameter. 

4. Sinusoidal variation of lattice parameter in the x 
direction gives subsidiary reflections at h ± a/L. In 
a cubic crystal there is nothing to distinguish the x direction 
from the y and z directions, and one would expect varia- 
tions in these directions also, either coexisting through the 
crystal, or taking place in separate regions. The investi- 
gations of Daniel and Lipson (1944) indicate that the latter 
is more likely. Points of the reciprocal lattice of the form 
//, k, I will thus be accompanied by three pairs of sub- 
sidiary reflections, the intensities being proportional to 
h 2 , k 2 , I 2 . Those of the form h, k, o will have only two 
pairs, since one of the indices is zero, and those of the form 
h, o, o will have only one pair. The separation of the 
subsidiary reflections from the main line in a powder 
photograph will not be a simple function of 0, but will 
depend on the indices of the reflection. This may be 
seen as follows. (Compare chapter VI, paragraph 2.) 
A reflection whose indices are h, k, I will reflect at an 
angle given by 

(2«/A) 2 sin 2 = h 2 + k 2 + I 2 . 

One whose indices are k + a/L, k, I will reflect at the 
slightly greater angle + 60 given by 

{2a/X) 2 sin 2 (0 + 60) = (h + a/L) 2 + k 2 + I 2 , 
or approximately 

(2a/A) 2 sin 6 cos 60 = ha/L, 
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The separation of the side-bands is thus greatest for reflec- 
tions of the form hoo, and least for those of the form hhJi. 
If h, k and / are all non-zero and unequal there will be 
three side-bands on each side of the main line, and their 
overlapping with the main line and one another makes 
measurements of intensities difficult. 

5. Equation (8) for the amplitude of reflection from 
a structure with sinusoidal variation of structure amplitude 
does not involve any approximations, but equation (6) for 
variation of lattice parameter depends on the assumption 
that Ldh'/a 2 is small compared with unity. An exact 
treatment hardly comes within the scope of this book, but 
it may be remarked that the second exponential in equation 
(4) is represented exactly by the series 

exp {i(L6h'/a*) cos {znna/L) } = J (Ldh'/a*) 
+ 21 J \(Ldh'/a 2 ) cos (zTtna/L) 

- iJ^Lbh far) cos (27m. 2a/ L) 

- 21 J z (L6h'/a 2 ) cos (27tn . 2a/ L) + ••; . (10) 

where J Q , J ly J 2 , '" are Bessel functions of the o ,h , i 8t 
2 nd , ••• order. Each reciprocal lattice point is therefore 
flanked by a series of side-bands at distances of ± a/L, 
± 2a I L, ± 3a/ L, •••, the intensities of the central reflection 
being proportional to J%, and the intensities of the others 

being proportional to J'j, J% ••'. It is a property of Bessel 
functions that 

JI + 2JI + 2JI + - = 1, . . (n) 

so that the total intensity is the same as that from an undis- 
torted lattice. (There is only one central maximum, but 
two side-bands of intensity J\, two of J'i, &c.) Daniel and 
Lipson (1944) find that the higher order side-bands do 
appear, and that when all factors are taken into account 
the intensities are in agreement with expectation. One of 
the difficulties is ' extinction ', a phenomenon that reduces 
the intensities of the main lines relative to those of the side- 
bands. The original paper must be consulted for full details. 
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6. X-ray diffraction from aluminium crystals containing 
4 per cent by weight (175 per cent by atoms) of copper has 
been investigated by Preston (1938) and Guinier (1938). 
The phenomenon is somewhat complicated ; at room 
temperature layers of copper segregate on planes perpen- 
dicular to the axes of the parent crystal, but at higher 
temperatures (about 200° C.) the copper layers ' evapor- 
ate ' and layers of Al 2 Cu form. The Al 2 Cu has not its 
usual structure, but one similar to that of fluorite (CaF 2 ), 
apparently because this structure can fit into the aluminium 
lattice with less strain. The sheets of Cu atoms at room 
temperature appear to be about 45 A square, 200 A apart, 
and about two atomic layers thick. The layers of Al 2 Cu 
produced at higher temperatures are considerably larger, 
and their diffraction effects sharper. No detailed treat- 
ment of the diffraction effects has been given in either case, 
but there is one very interesting phenomenon. The 
X-ray photographs show streamers from the ordinary 
reflections on one side only, in a way that indicates that 
the intensity of reflection in reciprocal space is not sym- 
metrical about the points of the reciprocal lattice. The 
h k I reflection has streaks running out towards 
h + i, k, I ; h, k + 1, / ; and h, k, I + 1 ; but not towards 
h — 1, k, I ; h, k - 1, / ; or h, k, I - 1. This is one of the 
few phenomena that are not effectively symmetrical about 
the points of the reciprocal lattice, and is at first sight 
rather puzzling. It is, however, qualitatively explained 
by a combination of varying lattice parameter and varying 
structure amplitude. Copper atoms are smaller than 
aluminium atoms, and contain more electrons, so that the 
scattering factor is high where the spacing is low, and vice 
versa. If both vary sinusoidally 



x n = na + (Ld/zna) cos (27tna/L), 
and F n = F[i +g sin (znna/L)] ; 



• (12) 



in order to ' synchronize ' the variations in the way postu- 
lated above the wave-length L must be the same in both 
8 
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x n and F n , and one must vary as the cosine, the other as the 
sine. The amplitude of reflection as a function of h' 
is then 

G(/0=^„exp{27n7^» . . (13) 

= FZ[\ + g sin (znna/L)] exp {2nih'[n + {Ld/zna-) 

n 

cos (27111a/ L)] }, 

which becomes, on expanding and neglecting products of 
small terms, 

G(h') = Flcxp {zTiinh'} 

+ iF[(L6h' /2a 2 ) -g/2]Zexp {2nin{h' + a/L)} 

+ iF[(Ldh'/2a*)+g/2]£lxp {z7cin{h' - a/L)} . . (14) 

The first of these sums gives a maximum of intensity 
proportional to F- at h' = h, where h is any integer, the 
second gives a maximum of intensity proportional to 
F 2 (L6h'/2a 2 -g/ 2 ) 2 for h' = h - a/L, and the third a 
maximum of intensity proportional to F 2 (Lbh' /2a 2 + g/2) 2 
for h' = h + a/L. It is clear that the inner maximum will 
be weaker than the outer, and will practically vanish if 
Lbh/a 2 is approximately equal to g. This simple model 
gives merely a pair of extra spots accompanying the main 
reflection, but it is not to be expected that L will be constant 
throughout the crystal, and a proper treatment would 
undoubtedly give a continuous streak extending from the 
main reflection towards h + 1, the intensity on the high- 
order side being greater than on the low-order side. 
Similar layers perpendicular to the jy and z axes will produce 
streaks running outwards towards k + 1 and / + 1. 

7. After suitable heat treatment (' annealing ') the spots 
or lines on X-ray photographs of metals are sharp, but if 
the metal is deformed by bending, twisting, rolling, draw- 
ing, hammering, &c, the spots spread out into Debye- 
Scherrcr lines, and the lines become broadened. Clearly 
the original perfection of the crystal lattice has been lost, 
but there has been some controversy over the nature of the 
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imperfection. There have been three main suggestions : 

(i) that the metal is broken up into ' crystallites ' so small 
(io~ 5 to io~ 8 cm. in linear dimensions) that dif- 
fraction broadening occurs ; 

(ii) that the metal is broken up into crystals (about io -4 
cm. in linear dimensions) with differing lattice 
parameters ; and 

(iii) that the crystals of the metal remain fairly large (about 
io -4 cm. in linear dimensions), but are elastically 
distorted. 

The distinction between the second and third suggestions 
lies chiefly in the nature of the distortion. According to 
the second the distortion is merely a uniform expansion 
or contraction of the whole crystal, whereas the third 
includes non-uniform expansion or contraction, as well as 
twists and bends. 

A careful study of the variation of the line-broadening 
with wave-length, Bragg angle, and indices of reflection 
can decide between the suggestions. The first would yield 
a broadening increasing with / sec 0, and possibly varying 
somewhat with the indices of reflection, if the crystallites 
tended to have a more or less regular non-spherical shape. 
The second would yield a broadening increasing with 
tan 0, but independent of A and the indices of reflection. 
The third would yield a broadening increasing with tan 6 
and independent of X, but it would be independent of the 
indices of reflection only if the tensile strain were, on the 
average, independent of crystallographic direction. Most 
crystals (tungsten seems to be an exception) are not elas- 
tically isotropic, so that strain independent of direction 
would imply a variation of stress with direction. On 
general considerations of static equilibrium it seems likely 
that stress, rather than strain, would tend on the average 
to be independent of direction, with a consequent variation 
of the strain, and a dependence of the line-broadening on 
the indices of reflection. 
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On these qualitative considerations the experimental 
evidence is unquestionably in favour of the third sug- 
gestion. It is found that the broadening increases with 
tan and not with sec (Brindley, 1940 ; Stokes, Pascoe 
and Lipson, 1943 ; Smith and Stickley, 1943), is inde- 
pendent of the wave-length of the X-rays (the same refer- 
ences, and Lipson and Stokes, 1943), and depends mark- 
edly on crystallographic direction (Wood, 1932 ; Stokes, 
Pascoe and Lipson ; Smith and Stickley). Though the 
experimental evidence makes it almost certain that lattice 
distortion is the main cause of the line-broadening, the 
theoretical problem of relating quantitatively the breadth 
of the lines and the properties of the crystal remains. It 
is as much a problem in the theory of elasticity as in X-ray 
diffraction, and possibly the nature of the glide processes 
by which excessive strains in the lattice are relieved must 
be taken into account also. No exact solution has been 
found, and perhaps each type of metal and each mode of 
deformation requires special treatment. It is, however, 
possible to give an approximate general theory. 

8. For calculating the line-broadening it is necessary 
to know Y(t), the mean value of FF* for two cells separated 
a distance t in the h k I direction (chapter VII, para- 
graph 8). In a distorted crystal the structure amplitude 
of a cell at (x, y, z) will differ from that at (o, o, o) for two 
reasons : the cell is displaced bodily, and it suffers a rota- 
tion. It is hardly possible to evaluate Y(t) without 
assuming a particular variety of distortion. It will, how- 
ever, be true in general that for large values of t the relative 
displacement of the cells will be large and random, and 
Y(t) will therefore vanish. For t small the displacement 
of one cell relative to another a distance t away will be 
e,„ the local value of the strain in the h k I direction, 
multiplied by t. The effect of the rotation on the phase 
of the scattered X-rays is small compared with that of the 
displacement, and may be neglected in a first approxima- 
tion. The structure amplitudes of the two cells will 
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therefore differ only through a phase factor 

exp {zniSe h t}, . . . (15) 

where S m 2 sin 6/X is the distance from the reciprocal 
lattice point in question to the origin.* The required 
function Y(t) is then the mean value of 

FF* = F 2 exp {27tiSe,,t}. . . (16) 

If the fraction of the crystal for which e h lies between e and 
e + de is y)f,(e) de, this mean value becomes 

Y(t) = F 2 [ "fa(e) exp {xntSet} de, . (17) 
J -» 

and the apparent particle size is 
e = [Y(o)]-i\ +a> Y(t)dt 

J —CO 

= 1 I i/» /; (c)exp {27tiSet}de\ dt 

= **T [P'toto exp {27ti(St)e} de^ d{St), 

= 5-%(o), (18) 

by the relations between Fourier transforms (chapter VII, 
paragraph 5). The integral breadth of the reflection is 

P = X/& cos =2 tan 0/0 A (o). . . (19) 

For a given reflection ft increases with tan 6 and is inde- 
pendent of X, as stated above, and 

A (o) = 2 tan 0/0. . . . (20) 

Measurement of integral breadths thus gives the fraction 

* In a cubic crystal the components of the relative displacement 
along the crystal axes are ei,t.h/\/N, eht.k/y/N, est.l/y/N, 
where N = h 2 + k s + I 2 . The difference in phase between the 
rays scattered by the two cells is therefore 
27r(h.e*th/VN + k.entk/VN + l.e„tl/VN)/a 

= 2Tr\/Neht/a = 277(2 sin 0/X)ei,t. 
The expression is, however, valid for non-cubic crystals (Stokes 
and Wilson, 1944). 
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of the crystal for which there is no strain in the hkl 
direction. In fact, it is fairly easy to show that when 
proper allowance is made for other causes of line-broad- 
ening the line profile is tft h ( - e cot 0) within a proportion- 
ality factor. The negative sign merely indicates that 
places where the lattice is expanded reflect at lower angles 
than the ideal, and places where the lattice is contracted 
at higher angles. 

Equation (20) may be regarded as a complete, though 
approximate, solution of the diffraction problem. There 
remains the problem of evaluating 1/^(0) in terms of the 
elastic, and perhaps the plastic, properties of the substance. 
One simple assumption is that 0,,(o) is inversely propor- 
tional to the maximum strain in the hkl direction, and that 
this is given by the maximum stress, assumed independent 
of direction, divided by the value of Young's modulus 
for the hkl direction. For a cubic substance this gives 

P cot = A + BH 

where A = 4/w s lv 

B = - 8p mux (s n - s u - 2s ti ), > ( 2I ) 

H = (h*k* + kT- + W)/<*" + k°- + /2)2 J 
/> inax . is the maximum stress and s n , s l2 , s u are elastic 
constants. Stokes, Pascoe and Lipson (1943) have found 
that this equation gives quite good agreement with experi- 
ment, and that p mAX is approximately the same as the 
known values of the yield stress (Megaw et al., 1944, 1945), 
but Smith and Stickley found only qualitative agreement. 
Stokes and Wilson (1944) have attempted to take shear 
stresses as well as direct stresses into account, using the 
approximation of paragraph 4, chapter VII, for the mean 
value of FF*. This gives 

Y(t) = <^exp {zniSe^t}} = F*exp {- 27i*SH\e}^} 
and s = [Y(o)]-^Y(t) dt = S^^))-*, 
P cot = 2 Vi^^f,})- 
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Thus ft cot is proportional to the root mean square value 
of the strain in the hkl direction. On the assumption that 
the stress is on the average spherically symmetrical this 
leads to 

(icotQ = ^/(A +BH) 

where A - S7i{(s n +2si.Xpn>-¥n( s i2+2s n Kp 2 12 >}\ , 22 \ 

b = - p*{(fa - *i 2 ) 2 - 4&K&> ; J 

<(/>'j,> and <(/>?•>) are the mean square values of the direct 
and shear stresses and H has the same value as before. 
The experimental results so far do not appear to be 
sufficiently accurate to decide whether (21) or (22) is the 
better approximation. 
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CHAPTER IX 

EFFECTS OF THERMAL MOTION 

i. The atoms forming a crystal are in a state of continual 
agitation about their mean positions. The period of vibra- 
tion of an atom is very small compared with the time 
required for taking an X-ray photograph, so that the 
photograph records the time average of the X-ray scatter- 
ing. The period of vibration of an atom is, on the other 
hand, large compared with the period of vibration of X-rays, 
so that at any instant the scattering is the same as if the 
atoms were at rest, but displaced from their mean posi- 
tions. For the interpretation of X-ray photographs it is 
therefore necessary to find the time average of the scatter- 
ing from a crystal in which the atomic displacements take 
on all values appropriate to its temperature. The calcu- 
lation falls naturally into two stages : (i) the determination 
of what functions of the atomic displacements are involved 
in the average scattering, and (ii) the determination of the 
average values of these functions in terms of the absolute 
temperature and the elastic constants of the crystal. 

2. Suppose that because of the thermal agitation the 
origin of the n ih unit cell is displaced from its mean position 
of .v„ = n x a, y n = nj), z n = n 3 c, to 

x' n = n x a + X, n I 



y' n = n 2 b + y n , > • • • CO 
z' n = n a c + Z n . J 

The effective structure factor of the cell will therefore be 
changed from F to 

F n - Fcxp {27ii{hX„/a + kYJb + IZJc)}. . (2) 

It is here assumed that the unit cell is displaced as a whole, 
which is not strictly true for crystals with more than one 
atom per unit cell, but it will be seen later that the assump- 
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tion is justifiable. The function determining the scatter- 
ing from any imperfect crystal is Y„„ the mean value over 
position in the crystal of F n Fn +m for two cells separated 
by the translation m x a, m 2 b, m z c (chapter VII, equation 
(24) ). Thus the problem is to find the mean value of 

F n F n+m = F'exp {27ti[h(X n - X n + m )/a + k(Y n - Y n+m )/b 

+ l(Z n - Z H+m )/c]}. . (3) 

The average value of this over all cells in a crystal at any 
instant is equal to its time average for a given pair of 
cells with the required separation. It was shown in chap- 
ter VII that the mean value of exp (ix) is approximately 
exp (- £<•**»» so tnat » 'f quantities belonging to the first 
cell are unprimed, and those belonging to the second cell 
are primed, the average value of equation (3) is 

Y m = F 2 exp {- 2n\[h(X - X')/a + k(Y - Y')/b 

+ i(Z - Z')A] 2 >}, . ( 4 ) 

where the angle brackets represent a time average. As a 
matter of fact the result of using this approximation is here 
correct, though the proof of its correctness is difficult 
(Born, 1942-3). The expression to be averaged is easily 
rearranged to 

(hX/a+kY/b+lZ/c) 2 +(hX'/a+kY'/b+lZ'/c)* \ ,s 

-2(hX/a+kY/b+lZ/c)(hX'/a+kY'/b+lZ'/c). J W 

The first two terms are always positive, being perfect 
squares, and will have the same mean value. Let 

M = 27i\(hX/a+kY/b+lZ/cy-), . (6) 

where M is necessarily positive. At first one might think 
that the average value of the third term would be zero, 
except for the special case m 1 = m 2 = m 3 = o, when it would 
cancel the first two. Cells close together in the crystal, 
however, will on the whole tend to move in the same way, 
the tendency diminishing with increasing separation of the 
cells. The average value of the third term will not vanish, 
therefore, but will start at - zM for the m's all zero, and 
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decrease numerically to zero with increase in the m's. Let 

P,n = Vt\{hX/a +kY/b+ IZ/c) 

{hX'/a + kY'/b + lZ'/c)>. . (7) 

Since this is symmetrical in the primed and unprimed 
letters P m = P_ m . Equation (4) becomes 

Y m -F*exp(P m -2M). . . (8) 
Since P„, decreases from zM to o as m increases, Y m 
decreases from its initial value of F 2 , approaching the con- 
stant value F 2 exp (- 2M) asymptotically. We have seen 
that a constant value of Y m means a sharp reflection, and 
a value decreasing to zero means a broadened reflection. 
In the present case Y m can be regarded as the sum of two 
parts ; a constant part F 2 exp ( - 2M), and a variable part, 
which it is convenient to write F 2 W m , decreasing from an 
initial value of .F 2 [i - exp (- 2M)] to zero with increasing 
m. Equation (8) becomes 

Y m =F 2 exp(-2M) + FW m> 
where W m - [exp (P,„) - 1] exp ( - 2M), 

and the average intensity of reflection from the crystal as 
a function of position in reciprocal space is (chapter VII, 
equation (24) ) 

H(u, v, w) m H x {u, v, w) + ll 2 (u, v, w) 
where 

!!,(,!, v, to) = IMF 2 exp (- zM)SV m 

m 

exp {- 27ii{tn{u + m 2 v + m a w)} 
and 
H£u t v t w)~U-*F*EV m W m 

m 

exp {- 27ii(m l u + m 2 v + m 3 w)}. 

The first of these is the ordinary expression for the inten- 
sity of reflection from a perfect crystal, except that it is 
diminished by the factor exp ( - 2M). The second is like 
that for the intensity of reflection from a crystal containing 
mistakes. The scattering from a crystal in thermal agita- 
tion consists therefore of two parts, the first the same, 
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(9) 



(.0) 



(») 



except for a reduction in intensity, as that from a perfect 
crystal with stationary atoms, and the second a broad re- 
flection similar to that from a crystal with mistakes. Since 
W m is equal to W_ m the maximum of the broadened 
reflection will coincide in reciprocal space with that of the 
sharp. The total intensity, the integral of the sum of 
the two parts over one cell of the reciprocal lattice, is the 
same as that from an undisturbed crystal. In general W m 
decreases very rapidly with m, and the broadened part of 
the reflection may be written 

// 2 («, v, w) = (V/Uabc)F*fflW(x,y, z) 

exp {- 27ii{xu/a + yv/b + zw/c)} dx dy dz . (12) 

where x = m^a, y = m 2 b, z = m 3 c and the integrations ex- 
tend from - 00 to +00. Within a factor, therefore, H 2 
is the Fourier transform of W. 

3. Logically the next step is to evaluate M and P m in 
terms of the absolute temperature and elastic constants of 
the crystal, but a digression on their physical significance 
and qualitative behaviour is perhaps desirable. The inter- 
pretation of M is the simpler, being 2n 2 S 2 times the mean 
square displacement of the cells in the hkl direction. This 
is easily seen for a cubic crystal, for then 

M = 2n\{hX/a + kY/b + IZ/c) 2 } 



= 271- 



Jt* + k 2 + / 2 



<(h'X + k'Y + l'Z) 2 } 



(13) 



= 27i 2 S\{h'X + k'Y +l'Z) 2 \ 
where h', k', /', the indices of reflection divided by 

are the direction cosines of the hkl direction. The expres- 
sion in parentheses is the component of the displacement 
X, Y, Z in the direction determined by h', k' } I', or, in 
other words, the displacement in the hkl direction. For 
crystals of other systems the geometry is more complicated, 
but the result is the same. 
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Since the intensity of the sharp component of the reflec- 
tion is reduced by the factor exp ( - zM), temperature 
motion reduces the intensity of all sharp reflections. The 
reduction is on the whole greater for the higher orders, 
that is, those with large values of S = 2 sin 0/X, but it is a 
smooth function of S only if the crystal is cubic or elastically 
isotropic. The mean-square displacement will be greatest 
in directions in which the restoring forces acting on the 
atoms are least, and the intensity of reflection from planes 
perpendicular to these directions will be lower than that 
from other planes. The effect is very marked in substances 
like graphite, in which displacements perpendicular to the 
layers are much easier than displacements in the layers. 

In a similar way it can be shown that P m is 4Jt 2 S 2 times 
the mean value of the product of the displacements in the 
hkl direction of two cells separated by the translation m x a, 
m 2 b, m z c. This will depend not only on the ease of dis- 
placement in the hkl direction, but also on the relation 
between this direction and the translation relating the cells. 
Thus P m and M will be large for the same points of the 
reciprocal lattice, but for a given point P.„, will be largest 
when the translation m^a, mj>, m 3 c lies in directions for 
which the crystal is not easily deformablc, since the cor- 
relation between the motions of the cells will be greatest 
for these directions, and it will be least when the trans- 
lation lies in directions for which the crystal is easily 
deformablc. The same is true for W M ; it will be large 
for the translation lying along relatively undeformable 
directions, and small for it lying along readily deformable 
directions. The intensity of the fuzzy part of the reflec- 
tion is related to W m by equation (12), whence it follows 
that the variation of H 2 (u, v, w) with u, v, w is the converse 
of that of W m with m v m 2) m^. Thus there will be streaks 
of comparatively strong intensity in directions for which 
the crystal is easily deformed, and regions where the scat- 
tered intensity is weak corresponding to directions for 
which the crystal is not easily deformed. This behaviour 



' 



is similar to that of the scattering from small crystals, dis- 
cussed in chapter III. The region of scattering is extended 
in directions for which the dimension of the crystal is least 
(face perpendiculars), and is contracted in directions for 
which the dimension of the crystal is greatest (corners). 
Small thickness corresponds to easy deformability and 
large thickness to difficult. 

4. The evaluation of M and W m bears many resemblances 
to the calculation of the heat capacity of crystals. The 
motion of the atoms can be analysed in terms of ' normal 
co-ordinates ', each of which is equivalent to a wave of 
definite frequency traversing the crystal. There are three 
times as many normal co-ordinates as there are atoms in 
the crystal, the number of ' degrees of freedom ' being the 
same whether the motion is regarded as made up of waves, 
or of motions of individual atoms. If there are more than 
one atom in the true unit cell, the waves divide into two 
groups, called ' acoustic ' and ' optical '. The former are 
of low frequency, with wave-lengths large compared with 
the cell dimensions, and the latter are of high frequency, 
with wave-lengths less than the cell dimensions. The 
number of waves in the acoustic group is 3N, where iV is 
the number of unit cells, the remaining waves being optical. 
The energy of the crystal is equal to the sum of the ener- 
gies possessed by the waves. It is known from quantum 
theory that the average energy of a vibration of frequency v 
at the absolute temperature T is 



E., = 



hv 



exp (hv/kT) - 1 
= kT for T»hv/k. 



+ Ihv 



■ (14) 



In this equation h is Planck's constant and k is Boltzmann's 
constant. If the number of waves per unit volume with 
frequencies between v and v + dv is Q,. dv, the total energy 
of a crystal of volume V is 

E-VfB&.dp, . . . (15) 
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the range of integration covering all frequencies. The heat 
capacity per unit volume is then C = V- l (dE/dT) r . The 
well-known Debye expressions for E and C„ follow from 
the assumption that the crystal behaves as an elastically 
isotropic continuum for frequencies up to a certain maxi- 
mum v„„ fixed by the requirement that the total number 
of frequencies should be three times the number of atoms. 
It is clear that this assumption applied to the problem of 
thermal scattering will obscure most of the interesting 
effects, since the mean-square displacement will be inde- 
pendent of direction, so that the weakening of the sharp 
maximum will depend only on S z and not on h, k, I, and 
the cloud of diffuse scattering around the sharp maximum 
will be spherically symmetrical and show no streamers or 
other peculiar effects. A detailed treatment of M and W m 
would, however, require an extensive knowledge of lattice 
dynamics, and we must confine ourselves to the Debye 
evaluation of M and a first approximation to the diffuse 
scattering. In order to keep the equations as simple as 
possible we shall assume that the axes arc orthogonal and 
that there is one atom per unit cell. 

5. A single wave of frequency v traversing the crystal 
will displace the n lh unit cell in the hkl direction by an 
amount varying sinusoidally with the time. Actually there 
are always three independent waves of a given wave-length 
— corresponding to the one longitudinal and the two trans- 
verse waves in an isotropic medium — and in the evaluation 
of M we may assume that one of these is responsible for 
the displacement in the hkl direction, and the other two 
are responsible for the displacements at right angles to it. 
The displacement of the n lh unit cell in the hkl direction 
is then of the form 



$,. = A v COS 27lVt, 



(16) 



where A,, is the amplitude of the wave of frequency v. 
To calculate M we want the mean value of £ 2 as a func- 
tion of the absolute temperature, and this is obtainable by 



way of the energy of the vibrating crystal. Since the mean 
value of cos 2 is one-half, <£*> = J<j4*>. The kinetic 
energy of the cell is \p$ = 27ch ,2 /.iA; sin 2 27ivt, where fi is 
the mass of the cell, so that the mean value of the kinetic 
energy is 7ih- z fi(A]^. The mean energy of the cell is twice 
this, the mean values of the potential and the kinetic 
energies being equal, so that the total energy of the crystal 
due to this wave is 

E„ = 2Nti 2 v*h(A% 
so that <| 2 > = \<j*> = Ej/4X*Nfa>*. . . (17) 

The mean-square displacement of a cell in the hkl direc- 
tion is given by equation (17) summed for all values of v 
that are effective, so that 

M = 27t 2 S 2 <£ 2 > 

= 27i*SKm*>VQ : .dv 

1 -^ dv, 

6N t u] v 2 

The factor of one-third comes in because only one-third 
of the waves are effective in producing displacements in 
the hkl direction. For an isotropic medium, as assumed 
in the Debye theory, Q„ = i2jrC-V, where C is an aver- 
age value of the speed of propagation of elastic waves. 
The derivation of this equation will be found in any good 
textbook on heat, for example, in chapter XX of Heat and 
Thermodynamics by J. K. Roberts. Thus on the assump- 
tions of the Debye theory 

where p is the density of the crystal. Since 



(18) 



M = 



(19) 



f 

Jo 



12.-TC- 3 \Y-dv = T,N, 



qnV 3 ifJiV 



C 3 =' + v* = 
3^ 



(20) 
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Substituting the value of E,. from equation (14) gives 



M = 



,f'Te 

.Jo L e 



2 /^mJo lexp (hv/kT) - 1 
or, with the substitutions = hv m /k, x 
3 k 2 S 2 T 



+ \hv \dv 
= 0/T,y=hv/kT, 



M = 



2/.tkQ°- 



where 



«A(*)=-f 

XJqCX 



#v) 

.y dy + x 
xp(:y)-i 4 






The integral occurring in i/»(.v) is not expressible in finite 
terms, but a table from which it may be obtained is given 
by Debye (1914) and in the Internationale Tabellen, p. 574. 
A selection of values of ifj(x) is given in table I. The 
quantity & is the same as the ' characteristic temperature ' 
occurring in the Debye theory of heat capacity. For x less 
than unity (temperatures greater than the characteristic 
temperature) ifj(x) is nearly constant, and 

-3M • • • « 

within three per cent. For greater accuracy in this range 
the expression 

#») = 1 + *'/36 • • • (23) 
may be used. 

6. Equation (11) for the diffuse intensity // 2 (m, v, w) is 
more or less rigorous, but it is difficult to proceed further 
with it without making approximations. At not too high 
temperatures exp (P,„) in equation (9) for W m may be 
replaced by 1 + P m , giving 

W m = e x P (- 2 M)P m . . . (24) 

One can then evaluate P m in terms of the amplitudes of 
the normal co-ordinates, much as was done for M, and 
obtain H 2 (u, v, w) from equation (12). The whole process 
really amounts to finding the Fourier transform of (A'*,} 
and then finding the Fourier transform of the result, giving 
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<^4J> again. The intensity at u, v, w is then proportional 
to <^4?.>, with v satisfying the relation 

v 2 /C 2 = « 2 /a 2 + v*/b* + w 2 /c 2 . 

The same result can, however, be obtained by a more 
elementary argument based on the treatment of a sinu- 
soidally-distorted crystal in the previous chapter, and as 

TABLE I 
Values of the function tfi(x) 



x 

mo 



o 
i-ooo 

6 
1-771 



1 
1-028 

7 
1-984 



2 
I-I07 



2-20S 



3 
1-233 

9 

2-433 



4 
1388 

10 
2-664 



S 
x-57i 



space does not permit including both methods the simpler 
will be given. 

From the point of view of the X-rays each elastic wave 
traversing the crystal is equivalent to a sinusoidal distortion 
of the lattice. It has been shown that a sinusoidal distor- 
tion in the x direction of amplitude Ld/2Jta and wave- 
length L produces a pair of satellite spots accompanying 
each reciprocal-lattice point at actual distances of ± l/L 
from it and of intensity (Ldh/za 2 ) 2 smaller. This result 
has now to be generalized in two ways. If the wave normal 
is not parallel to the x axis, the satellite spots will be dis- 
placed from the point of the reciprocal lattice along a line 
parallel to the wave normal, not along a line parallel to 
the x axis. Secondly, it is the component of the displace- 
ment in the kkl direction that produces the scattering round 
the point hkl of the reciprocal lattice, so that the effective 
amplitude of the wave is its actual amplitude multiplied 
by the cosine of the angle between the hkl direction and 
the displacement direction. A wave of amplitude A and 
length L will therefore produce satellite spots of intensity 
proportional to (nFAS) 2 , where r is the cosine of the angle 
between the displacement direction and the hkl direction, 
and S = 2 sin 0/?. is the distance of the point hkl from the 
9 
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origin of the reciprocal lattice, the spots lying at distances 
of ± i/L, measured along the wave normal, from the point 
hkl. This reduces to the earlier result for a wave of 
amplitude Ld/27ta with its wave normal along the x axis. 
For a given value of L there are three waves of different 
frequencies, speeds, amplitudes, and directions of displace- 
ment. The intensity at a distance i in a given direction 
from the point hkl will thus be proportional to 

7i*s*(r*Ai+riAi + riAi), . . (25) 

where quantities belonging to the three waves are distin- 
guished by the subscripts i, 2, 3. In calculating M we 
could assume that one of the three waves produced a dis- 
placement in the hkl direction and the other two at right 
angles to this direction, but here we have not assumed 
isotropy of the crystal, and all three must be taken into 
account. The average values of the squares of the A's 
are given by equation (17), so that the intensity is pro- 
portional to 

2N{i\_ v\ v\ v\ J 

s* vr\E x m 2 m~-\ 
2N{is*l c\ ci c 2 y 

where the v's and C's are the frequencies and speeds of the 
three vibrations of wave-length L. This expression gives 
the intensity at s as a fraction of the intensity at the point 
of the reciprocal lattice, which is N 2 F 2 , or N 2 F 2 exp ( - zM) 
when allowance is made for the temperature motion. The 
three waves travelling in the opposite direction will con- 
tribute an equal amount, so one would expect the intensity 
at u, v, w to be given by 
Ho(u, v, to) 




NF 2 exp (- zM) S* [7^ /% £]£,! 
V'LCf Ci C| J* 



(27) 



where s 2 = u 2 /a 2 + v 2 /b 2 + w 2 /c 2 for an orthogonal lattice. 
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The correctness of the proportionality factor can be veri- 
fied by showing that the total diffuse intensity is equal, to 
the first order in M, to the intensity lost by the central 
spot. On the Debye theory there is a maximum value of v, 
independent of direction in the crystal, and hence the cloud 
of diffuse scattering is spherically symmetrical, and stops 
abruptly at a limiting value of s. For the purpose of check- 
ing the proportionality factor we may ignore the differences 
between the E's and C's, and take the total diffuse intensity 
to be the integral of equation (27) over the interior of a 
sphere of reciprocal space of radius v m /C. The sum of 
the J"s will be unity under these conditions, and the total 
diffuse intensity becomes 

JJJ// 2 (m, v, w) du dv dw 

= abcMH 2 {u, v, zo)d{u/a)d(v/b)d(w/c) 

= u[ mC H 2 {s).^is 2 ds 
Jo 

NUF 2 exp (-2M) fWO B, 

H ' Jo cv 



$ns 2 ds, 



or, since s =• i/L 



v/C and NU = V, 
'C 3 "Jo 



VF 2 exp(- 2 M) S2 



E„dv 



= F 2 exp{-2M).zMN, . . (28) 

by equation (19). The reduction of the total intensity of 
the central spot is 

7VT 2 -iVT 2 exp(-2M), . . (29) 

which is the same as equation (28) to the first order in M. 
Agreement to higher orders cannot be expected, because 
of the approximations made in deriving equation (27). 

7. Equation (27) shows that the intensity of the thermal 
scattering increases with increasing order of reflection (iS 2 ), 
decreases with increasing distance from the points of the 
reciprocal lattice (1 /s 2 ), and for temperatures greater that 
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varies like S-Texp (- KS 2 T), where Kis independent of T 
and (to the Debyc approximation) the indices of reflection. 
It would thus have a maximum value for S Z T = i/K. The 
approximations made, which are equivalent to the expan- 
sion of exp (/',„) to obtain equation (24), would hardly be 
valid for such large values of S 2 T, but it seems likely that 
the predicted behaviour is qualitatively correct, and that 
the intensity of reflection in the regions near the points 
of the reciprocal lattice does ultimately decrease with in- 
creasing temperature or order of reflection (Born, 1942-3, 
p. 317). The factor in brackets, involving the speeds of 
propagation of the vibrations, is in general strongly depen- 
dent on the direction of 5 and on the point of the reciprocal 
lattice concerned. It has been worked out for cubic 
crystals by Jahn (1942), but the expressions are too com- 
plicated to be quoted and discussed here. 
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